_:   i  sr«»trfhcei,tkwy9A3,H.y 
NEW  YORK  UNIVERSITY 
WASHINGTON  SQUARE  COLLEGE 
MATHEMATICS  RESEARCH  GROUP 

RESEARCH  REPORT  No.   172-4 


Propagation  of  Dipole  Radiation 
Through  Plane  Parallel  Layers 

by 
JEROME  LURYE 


This  Report  was  Prepared  and  Written  under 

Contract  No.  W28-099-ac-172 

tvith 

Army  Air  Forces,  Air  Materiel  Command 

Watson  Laboratories 

Red  Bank,  N.  J. 

November,  1947 


Hew  York  University 

Washington  Square  College 

Hatheaatics  Research  Group 

Research  Report  I72-U 

under 

Contract  l^o.   W2S~099-ac-172 


PROPAGATION  OF  DIP0L3  Rft.DlAiIOi:  THROUGH 
PLAilS  PARA.LLSXJ  LAYERS 

by 
Jerome  Lurye 


Hollis  R.    Cooley 
Technical  Sditor 


Title  page 
\  .,.  50  numbered  pages 

ri  }/''>  ilovemher  I9U7 


Morris  Kline 
Project  Director 


coi:TSi?rs 

Article  Page 

Abstract  2 

1.  Introduction  3 

Section  I.  lormulation  of  the  Protlem. 

2.  Description  of  the  Problem  7 

3.  Preliminary  Formulation  as  a  Scalar  Problem  g 
U.  Introduction  of  Cylindrical  Coordinates  10 

5.  Precise  Formulation  of  the  Problem  in  Terms  of  the  Hertz  Potential   11 

Section  II.  Dipole  outside  the  I'iultilayer 

6.  Description  of  the  Plane  Wave  Problem  lU 

7.  Mathematical  Expression  for  Plane  Waves  I7 
g.  Solution  of  the  Plane  Wave  Problem  I7 
9.  Summary  of  the  Plane  Wave  Problem  22 

10.  The  Superposition  Principle  23 

11.  Solution  of  the  Propagation  Problem  \/ith  the  Dipole  outside 

the  I'iultilayer  28 

12.  Validity  of  the  Generalized  Superposition  Principle  32 

13.  Conclusion  to  Section  II  33 

Section  III.  Dipole  within  the  Multilayer 

lU.  Statement  of  the  Problem  3U 

15.  Decomposition  into  Subsidiary  Problems  35 

16.  Solution  of  the  First  Subsidiary  Problem  35 

17.  Solution  of  the  Second  Subsidiar,,  Problem  39 

18.  Conclusion  of  Section  III  111 

Section  IV.   A  Special  Case 

19.  Statement  of  the  Problem  L1.3 

20.  Solution  of  the  First  Subsidiary  Problem  kk 

21.  Solution  of  the  Second  Subsidiary  Problem  U6 


2. 


Abstract 

This  paper  is  concerned  v/ith  the  problem  of  finding  the  strength  of  the 
electromagnetic  field  at  any  distance  from  a  dipole  source  under  certain  special 
conditions,  of  v;hich  the  principal  one  is  that  the  transmitting  medium  is  divided 
into  horizontal  layers  above  a  plane  earth,  uith  6,  /x,  and  (5"  constant  in  each 
layer.   Since  the  source  is  a  dipole  it  is  possible  to  express  the  field  in  terms 
of  a  Hertz  potential  function  with  one  component  only.  The  final  statement  for  the 
field  strength  is  in  terms  of  complex  integrals. 

There  is  first  considered  the  case  of  a  multilayer  contained  between  two 
serai-infinite  homogeneous  media,  A  stream  of  plane  \/aves  of  constant  amplitude  and 
polarized  so  that  the  H-vector  is  parallel  to  the  xy-plane  impinges  on  the  multilayer 
from  below,  and  the  field  in  the  inultilayer  is  obtained.   In  the  next  stage  of  the 
problem  a  dipole  source  of  radiation,  considered  as  a  superposition  of  plane  v/aves, 
is  assiuned  to  be  located  in  the  homogeneous  region  belov;  the  multilayer,  and  the 
field  due  to  the  dipole  is  obtained.   In  the  final  stage  of  the  problem  the  dipole 
is  located  within  the  multilayer  and  the  field  is  again  found.   As  an  illustration 
of  the  theory  a  special  case  is  treated  in  detail.   This  case  consists  of  a  layer 
of  finite  height  containing  the  dipole,  located  betv/een  the  earth  and  a  semi-infinite 
layer  above. 

The  complex  integral  form  of  the  final  result  may  permit  api:»roximations 
in  the  form  of  asymptotic  expressions  v/hich  \.'ould  have  computational  advantages  over 
the  conventional  forms  involving  slowly  converging  series.   ;/hile  this  theory  depends 
upon  the  assvmrption  of  a  flat  earth,  the  curvature  of  the  earth  is  taken  into  account 
by  means  of  the  modified  index  of  refraction. 


3. 

1.  Introduction 

The  object  of  this  report  is  to  explain  a  rather  general  method  for  deal- 
ing with  a  certain  class  of  radio  wave  propagation  problems  pertaining  to  non- 
homogeneous  atmospheres.  This  method,  which  consists  briefly  in  the  decomposition 
of  the  antenna  radiation  into  plaine  waves  and  the  subsequent  solution  of  a  plane 

v;ave  oropagation  problem,  is  by  no  means  essentially  new  -  it  has  been  used  by 

(1)*         (2)(^)(U) 
Weyl    and  others    ^    to  deal  \,'ith  special  cases  of  the  propagation  problem  - 

but  in  this  paper  we   shall  be  concerned  with  its  application  to  the  solution  of  a 

more  general  kind  of  problem.  The  method  possesses  several  features  worth  noting: 

a)  The  class  of  problems  with  which  it  deals  is  large  and  includes  all 
those  in  which  the  atmosphere  is  horizontally  stratified  in  layers  of  finite  thick- 
ness. - 

b)  The  solution  is  obtained  in  the  form  of  one  or  more  complex  integrals 
which  must  then  be  evaluated  by  some  suitable  procedure  ,-such  as  series  representa- 
tion»  asymptotic  expansion,  or  reduction  to  a  knovm  form.   This  situation  is  to  be 
contrasted  uith  the  one  v/hich  exists  when  the  more  usual  approstch  to  the  propaga- 
tion problem  is  employed,  namely,  that  of  expansion  in  so-called  "normal  modes." 
(For  a  comprehensive  description  of  this  procedure,  see  a  Radiation  Laboratory 
report  by  Furry   .  Also  papers  by  Watson  "^  ,  van  der  Fol  and  Bremmer   , 
Sckersley^' '^   ,  Pekeris^  \   Pearcey    ,  and  Freehafer^   .  )  These  normal  modes 
are  actually  transcendental  functions,  an  infinite  series  of  v/hich  is  used  to  repre- 
sent the  solution  to  the  propagation  problem.   Since  this  series  is  generally  very 
slowly  convergent,  the  effort  and  labpr  involved  in  obtaining;  a  numerical  result 
are  considerable.  It  is  then  natural  to  ask  whether  the  complex  integral  form  of 
solution  in  which  our  method  eventuates  is  any  easier  to  interpret  numerically  than 
is  the  normal  mode  form.  While  for  reasons  given  below,  we  have  not  in  this  paper 
investigated  the  practical  evaliiation  of  our  results  -  our  solution  is  therefore  a 
purely  formal  one  thus  far  -  it  seems  probable  that  corresponding  to  certain  physi- 
cal conditions  such  as  particular  kinds  of  variation  of  refractive  i»dex  v/ith  height, 
particular  ranges  of  the  ratio  of  wavelength  to  distance  from  the  antenna,  etc., 


References  cited  in  the  introduction  are  listed  on  page  6  at  the  end  of  article  1. 


our  complex  integral  expressions  would  require  louch  less  labor  for  their  interpreta- 
tion than  v/ould  the  corresponding  normal  mode  expressions. 

c)  In  this  paper,  our  method  is  used  in  conjunction  v;ith  a  single  compon- 
ent or  scalar  Hertz  potential  formulation  of  the  problem,  i.e.,  the  unknovm  will  be 
a  single  scalar  function  of  position  and  time,  just  as  it  is  in  the  normal  mode  pro- 
cedure.  Ue  may  also  make  a  comparison  in  this  regard  between  the  present  paper  and 
one  v;hich  is  to  appear  subsequently  by  R.  K.  Luneberg  of  this  project,  in  i/hich  two 
unknown  scalar  functions  (not  the  Hertz  potentials,  ho\/ever,)  are  employed.  The 
use  of  tv/o  scalar  potential  functions  broadens  the  class  of  problems  which  can  be 
treated  in  that  it  permits  one  to  consider  problems  i;ith  arbitrary  antennas;  more- 
over, Luneberg  defines  these  functions  so  that  they  are  applicable  to  the  situation 
in  which  the  atmospheric  refractive  index  varies  continuously  \;ith  the  height  in- 
stead of  stepv/ise  as  in  the  present  paper;   Luneberg' s  v;ork  is  therefore  more  gen- 
eral than  the  present  paper  in  the  two  respects  mentioned.  On  the  other  hand,  it 
should  be  pointed  out  that  the  present  paper  could  be  generalized  through  the  intro- 
duction of  a  second  Hertz  potential  function  to  include  antennas  as  arbitrary  as 
those  v/ith  v/hich  Luneberg 's  treatment  is  capable  of  dealing;   the  extension  of  our 
method  to  the  case  of  continuously  variable  refractive  index  is  not  so  immediate, 
hov/ever. 

Because  of  their  greater  generality,  Luneberg's  results  will  in  all  like- 
lihood be  investigated  with  respect  to  practical  interpretation  before  those  of  the 
present  paper.  It  is  for  this  reason  that  we  have  not  as  yet  undertaken  any  ajialysis 
of  the  complex  integrals  obtained  here,  preferring  to  wait  until  the  class  of  prob- 
lems for  v;hich  Luneberg's  procedure  gives  practical  results  has  been  decided  upon. 
The  present  report,  therefore,  is  in  no  sense  a  final  solution  to  the  problems  v/ith 
which  it  deals,  but  rather  constitutes  one  theoretical  approach  to  such  problems. 

In  this  paper,  as  in  all  recent  investigations  of  the  problem  of  propaga- 
tion of  high  frequency  waves  in  non-homogeneous  atmospheres,  v/e  shall  suppose  the 
surface  of  the  earth  to  be  flat  and  correct  for  this  departure  from  the  true  situa- 
tion by  a  now  standard  device  known  as  the  modified  index  of  refraction.  It  is 
commonly  assumed  that  the  atmosphere  is  stratified,  that  is,  that  the  index  of  re- 
fraction varies  only  with  hei^t  above  the  surface  of  the  earth.   ¥e  then  replace 
the  tme  variation  of  the  index  of  refraction  as  a  function  of  the  height  above  the 
surface  of  the  earth  by  a  uiodified  index  of  refraction  v/hich  is  given  by  the  formula 

^'^"^^     bn(b)    • 
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v/here  n(r)  represents  the  actual  index,  r  is  the  distance  from  the  earth's  center, 
and  b  is  aay   fixed  value  of  r,  usually  taken  to  be  the  earth's  radius  a.   If  \;e 
no\/  replace  r  by  a  +  z  this  formula  gives  the  modified  index  of  refraction  as  a 
function  of  z  alone,  z  being  the  height  above  the  earth's  surface  in  a  rectangular 
coordinate  system.   In  this  coordinate  system  x  and  y  no\i   represent  distances  along 
the  flat  Surface,  and  any  point  on  this  surface  corresponds  to  the  point  on  the 
earth  having  the  same  distance  and  direction  from  the  origin  along  a  great  circle 
route. 

Physically  expressed,  the  justification  for  the  modified  index  of  re- 
fraction is  roughly  the  follov/ing;  that  instead  of  having  the  earth  curve  away  from 
the  radiation  the  index  is  changed  in  such  a  v/ay  as  to  make  the  rays  curve  away 
from  the  earth.  The  mathematical  justification  for  the  use  of  this  modified  index 
of  refraction,  uhich  corrects  approximately  for  the  flattening  of  the  earth,  is 
given  most  clearly  in  a  paper  by  Pekeris   -^  ,  v/herein  estimates  of  the  validity  of 
this  approximation  are  also  given.  Essentially  the  approximation  is  valid  out  to 
ranges  of  the  order  of  the  radius  of  the  earth  and  to  heights  of  the  order  of 
1,000  feet  for  uave  lengths  larger  than  a  centimeter.   If  the  range  or  height  is 
exceeded  or  the  v/ave  length  is  decreased  the  approximation  becomes  poorer.  Actual 

data  are  given  by  Pekeris  in  this  report.  The  earth-flattening  approximation  is  con- 

(lli) 
sidered  for  purposes  of  ray-tracing  by  Freehafer. 

It  is  also  necessary  in  the  case  where  the  earth  is  treated  as  an  imper- 
fect conductor  to  correct  for  the  index  of  refraction  in  the  earth.   This  means 

that  in  place  of  the  relation  n  =  JEm     v/here  8.  and  Mare  the  dielectrit:  constant 

and  penaeatility,  which  holds  for  dielectrics,  one  must  use  n  =  >/^^  ij  v/here 

c     /     i  C 

t.  -  c  -  -r-r-  .   In  case  of  imperfect  condiictivity,  therefore,  it  is  necessary  to 


r.2  2 

r      r    ,       r 


replace  £  and  C  by  — 5 — r^- £.  and  — :r— -— -  (J  respectively,  so  that  the 

b^  [n(b)/  b^  [n(b)]2 

modified  index  of  refraction  results.   So  correction  is  made  for/U  since  it  is  unity 

in  practically  all  atmospheres  near  the  surface  of  the  earth  and  for  short  distances 

in  the  earth.   In  the  theory  of  this  paper,  v/hich  presupposes  a  flat  earth,  f .  and  C 

have  these  altered  values  wherever  they  appear. 
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Section  I 
Formulation  of  the  Propagation  Problem  to  be  Solved, 


2.   Description  of  the  Problem. 

We  shall  nov/  describe  the  kind  of  problem  to  which  our  method  is  applic- 
able. Referring  to  Fig.  1,  the  line  G  represents  the  surface  of  the  earth,  v/hich 
we   v/ill  assujne  throughout  this  report  to  be  perfectly  plane.  The  region  below  G, 


^  Ht/  >y^^i^*tJ  ^hr* 


^>J^'><^> 


^.i/^o<^0  Fig.  1 

YiOte;^  ,AA,<r   represent  the  dielectric  constant,  magnetic  permeability,  and  electri- 
cal conductivity,  respectively. 


of  course,  cons  tit  ites  the  earth  itself,  v/hich  will  generally  be  assumed  to  have  a 
finite  electrical  conductivity.   AboVe  G,  the  atmosphere  is  taken  to  be  horizontally 

stratified  in  layers  of  finite  thickness   d^ ,  d^,  ....  d  as  shovm.   The  boundary 

1   c:        n 

surfaces  of  each  layer  are  horizontal  planes,  and  \/ithin  any  one  layer  the  atmosphere 
is  perfectly  homogeneous,  undergoing,  however,  an  abrupt  change  in  its  electrical 


properties,  e.g.,  its  dielectric  constant  and  electrical  conductivity  from  one  layer 
to  the  next.   The  number  of  layers  n  is  finite,  and  after  the  n   layer  the  atmos- 
phere is  permanently  homogeneous  as  ue  move  upv/ard.   The  antenna  A  is  assumed  to  he 
a  vertical  electric  point  dipole  and  may  he  located  anyv/here  in  any  of  the  layers. 
The  problem  is  to  find  the  values  of  the  electromagnetic  field  intensity  vectors, 
E  and  K,  at  any  point  above  or  on  the  earth.  Before  proceeding  to  the  solution  of 
this  problem,  however,  \re   mast  first  formulate  it  somev/hat  more  simply  and  pre- 
cisely. 

As  stated  above,  this  pjropagrtion  problem  is  a  vector  problem,  i.e.,  we 
seek  the  field  vectors  3  and  H  as  functions  of  position  (the  time  variation  is 
ass\Maed  harmonic  throughout),  and  since  the  vector  unknowns  are  two  in  number,  in 
general  we  have  to  find  the  six  scalar  fixnctions  S  ,  E  ,  £  and  H  ,  H  ,  H  ,  i.e., 
the  six  components  of  E  and  H  in  tlie  directions  of  the  x,y,z-axes  shown  in  Fig.  1 
(the  y-axis  points  into  the  paper.)  If  now  \/e  could  find  some  property  of  our 
propagation  problem  \rhich  \;ould  enable  us  to  reformulate  it  so  that,  instead  of 
having  to  find  six  scalar  functions  v;e  need  find  only  one  in  order  to  solve  our 
problem,  \/e  would  obviously  effect  an  iraraense  simi^lif ication.   We  desire,  in  other 
v/ords,  a  single  scalar  potential,  from  v;hich  both  of  the  vector  solutions,  S  and  H, 
can  be  derived.   For  the  kind  of  problem  that  is  represented  by  Fig,  1  such  a 
scalar  potential  exists,  and  in  \/hat  follows  \/e  indicate  how  the  problem  is  formu- 
lated in  terms  of  it  and  ho\7  the  field  vectors  are  derived  from  it. 

3.  Preliminary  Formulation  as  a  Scalar  Problem. 

_♦    — * 
First  of  all,  let  us  consider  the  field  distribution,  S  and  H  ,  to  which 

the  antenna  A  v;ould  give  rise  if  it  v/ere  located  alone  in  a  homogeneous  medium  of 

infinite  extent  in  all  directions  and  characterized  by  the  electromagnetic  constants 

£  y    u.y    (J.   of  the  j   layer  of  Fig.  1,  this  layer  being  the  one  in  which  the 

antenna  of  our  general  propagation  problem  is  assumed  to  be  situated.   For  an 
electric  dipole  in  such  an  infinite  homogeneous  medium,  it  is  known  -^  that  the  H 
vector  has  no  component  in  the  antenna  direction,  i.e.,  with  the  orientation  shovm 
in  Fig.  1,  H  =  0  everywhere.   Moreover  -  and  this  is  the  important  point  -  if  the 
earth  and  the  stratified  atmosphere  are  nov/  restored  so  that  the  "homogeneous 


(15) 

Stratton,  J.  A.  :  Electromagnetic  Theory,  pp.  ^3^,1+35. 
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-*   — *  _    _ 

medium  field  vectors"  E  and  H  are  changed  to  the  vector  distriMtions  3  and  H 

v/hich  constitute  the  final  solution  to  our  propagation  problem,  then  this  solution 

field  is  also  characterized  by  the  condition  H  =  0.  In  short,  the  solution  to  the 

propagation  problem  illustrated  in  Fig.  1  is  a  so-called  "transverse  magnetic" 

field  \fith  respect  to  the  z  or  antenna  direction. 

The  latter  statement  follov/s  from  the  fact  that  the  boundary  surfaces 
of  the  various  media  are  all  horizontal  planes  perpendicular  to  the  antenna  ajcis 
and  hence  parallel  to  the  vector  H  .  Now  the  Hajn/ell  boundary  conditions,  which 
the  solution  vectors  S  and  H  must  satisfy  at  the  aforementioned  surfaces,  require 
only  tliat  the  component  of  E  parallel  to  the  surfaces  and  the  component  of  MH 
perpendicular  to  the  surfaces  vary  in  a  continuous  manner  from  one  side  of  each 
surface  to  the  other;   therefore,  since^K  is  already  parallel  to  these  surfaces, 
i.e.,  possesses  no  perpendicular  component,  it  is  evident  that  the  boundary  condi- 
tion on  ^H  may  be  met  by  assuming  that  the  perpendicular  component  of  yUH  is  equal 
to  zero,  and  hence  also  that  of  H.  Thus  the  assumption  that  the  solution  field  is 
transverse  magnetic  is  consistent  with  the  boundary  conditions  which  it  must  satis- 
fy. 

This  transverse  oagnetic  property  of  the  solution  to  o\ir  propagation 

problem  enables  us  to  formulate  the  problem  in  terms  of  a  single  unknown  scalar 

(l6 ) 
potential  function.  V/e  simply  make  use  of  the  v;ell  knovm  fact    that  if  in  a 

homogeneous  medium  characterized  by  electromagnetic  constants  £  ,yt<  ,  (5",  tliere  is 

an  electromagnetic  field  which  is  transverse  inagne^ic  \/ith  respect  to  the  z-direc- 

tion  of  some  xyz-coordinate  system,  then  a  scalar  function  of  position  and  time 

V|/(x,y,z,t)  exists,  having  the  following  two  properties:        "^  "' 

a)  W  satisfies  the  wave  equation 

V=Y-^e^-;..^=  0         (X, 

b)  The  components  of  the  E  and  H  vectors  of  the  field  are 
derived  from  ^  by  the  formulas 

X  <5zox'  y  dzay'  z  ^.2  ,2^' 

c>x  ^y 


(2) 


\-  ^-h  ^l^^-^l^".  \--'4^^%^  -<^W-  ^3) 


ITU) 


H,  =  0 


Stratton,    op.    cit.,   pp.    3IJ.9,    ff. 


■.;,  u 
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and  from  expressions  (2)  and  (3),  ^  and  H  satisfy  Maxwell's  equations  automatically 
"because  of  equation  (l). 

RemenlDering  that  v/e  are  considering  only  those  fields  which  are  harmouio 
in  tine  and  that  v/e  may  therefore  write  H  =  ^e  where  4*  -   H'(x,y,z),  a  func- 

tion of  position  alone,  we  see  that  a)  and  t )  "become:   , 

a')  ^^^J  +  k^^>=  0,  (u) 

where  k^  =    W  e  cj^  +  iu  a  CO     • 

^     azax        ^     azay  3x  3x'^ 

H     =-isi-     ^     .        H     = ^     J^     ,      H       =   0.  (6) 

Here  E  t  ^  »  etc.  now  represent  only  the  spatially  varying  parts  of  the 
X   y 

field  components,  the  field  components  including  the  time  variation  "being  given  "by 

-itOt   _^   -iUJt 

JL     e  ,  i  e      ,  etc. 

X       •  y        • 

U-.  Introduction  of  Cylindricnl  Coordinates. 

Our  propagation  problem  will  "be  somewhat  simplified  "by  the  use  of 
cylindrical  coordinates   v ,  if  ,   z,  where 


=  Jl^  -   "2 


'X  +  y 


(P=   arc  cos  ^  ,  (7) 

z  =  z  . 

If  the  antenna  of  ?ig.  1  is  located  on  the  z-axis  of  our  coordinate 
system,  then  it  is  evident  from  symmetry  that  the  resultant  field  is  independent 
of  the  angle  (4)  .      For  purposes  of  simplicity,  therefore,  v;e  shall  altvays  take  the 
antenna  to  be  in  this  position.   Kence  we  need  concern  ourselves  in  this  section 
only  with  potential  functions  of  the  form  vjiCr,  z),  and  the  components  of  E  and  H 
with  respect  to  the  cylindrical  coordinate  directions  are  given  "by 


* 

''Stratton,  op.cit.  p.  35O,  forms  (3)  3,   (U) 
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2,2 


H  =  0,     H.rt   =  -  -T-^  ^^   ,    H  =  0  (9) 

ijhich  are  analogous  to  equations  5  a^^^L  6. 

The  function  ^  is  knovm  as  the  Hertz  potential  of  the  electromagnetic 
field  in  question.  Let  us  nou  fornulate  the  propagation  problem  of  Fig.  1  in  terms 
of  the  Hertz  potential. 

5 .  Precise  Formulation  of  the  Pro"blem  in  Terms  of  the  Hertz  Potential. 

The  mathematically  precise  statement  of  the  propagation  problem  is  that 
ue  seek  a  function  VL/(r,z)  uhich  satisfies  the  follov/ing  conditions; 

A  -  Denoting  the  values  \/hich  the  function  (^   takes  on  in  the  j   layer 
of  Fig.  1  by  H^  . ,  then  as  long  as  the  antenna  is  not  located  in  the  j    layer,  we 
have  from  (U), 

P  P 

V   ^C  .  +  k.  ^J.  =0  every\;here  in  the  layer,      (10) 
where  kj  = /.  .  £.  (j^     -  i  ^  .   <S  .  u  . 

B  -  If  the  antenna  is  situated  in  the  j   layer,  then  l^.  again  satisfies 
equation  (10 )  everyv/here  in  the  layer  except  at  the  point  viiere  the  antenna  is 
located.  At  this  point  ^  becomes  infinite  for  the  following  reason.   In  the  immedi- 
ate vicinity  of  the  antenna  the  field  must  behave  like  the  "homogeneous  medium 
field"  of  the  antenna,  described  in  article  2.  For  an  electric  dipole  located  in 
an  infinite  homogeneous  medium  with  constants  £. . ,  u  .,   <5~ • ,  it  can  be  proved  that 

-*  -*  (18)  t  T 

the  field  E,  H  ,  can  be  characterized  by  a  Hertz  potential  having  the  form 

■   *  _   e^^*'  ^ 

^  j  ~     "^1  (11) 


'    This  condition  on  the  field  evidently  follo\/s   from  physical  considerations. 
IS.Stratton,    op.    cit.,   pp.    U3I-U33. 

T  T  We  are  usirig  M.K.S.   units   and  ss  a  consequence,  Vi>  .   should  really  be  multiplied  by 

1  "^ 

a  factor      >  .      For   convenience,    however,    \ie  liave   suppressed  this   factor   in 

J 
(11)   and  throughout   the   rest    of   the  paper. 
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where  R  is  the  distance  from  the  otservation  point  to  the  antenna.   Our  second  t^n- 
dition  on  UJ  . ,  therefore,  may  be  written 

1  ici        _  V  f  M  n  \ 

R  ->  0  ^'^     *^  , 

and  from   (11)   we  see  that    \i;.   must  tecome   infinite  at    the  antenna   in  consequence. 

C   -  Our   third  condition  on\s|/has   to   do  with   its   hehavior  at   large   dis- 
tances from  the  antenna.      We  require   on  physical  grounds  that  as   we  proceed  further 
and  further  from  the  antenna,   the  field  should  tend  to   zero,    i.e., 

^^''     KJ=  0.  (13) 

H  ->   00 

Furthermore,   we  require   that   at   sufficient   distances    from  the   dipole, 
the   field  should  "behave   like  an   "outgoing  v/ave",    i.e.,    a  \/ave   traveling  avmy  from 
the  dipole.      From  a  physical  point   of  viei;,    the  necessity  for  this   condition  is 
apparent  since  its  violation  \;ould  imply  the  existence   of  sources    of  energy  other 
than  the  antenna  itself,    which  is  contrary  to   our  assumptions.      The  mathematical 

statement   of  this   outgoing  wave  requirement,  knov/n  as   Sommerf eldfe   "radiation  condi- 

(19) 
tion",    is  that  V|/ should  satisfy  the   following: 


T?    - 


ii°     R(aJ^    _  iku;)  =  0.  •  (11+) 

>    00 


D  -  Our  final  conditions  on  ki/ relate  to  its  "beliavior  on  the  houndary  sur- 
faces of  the  various  strata.   This  "behavior  is  determined  from  the  electromagnetic 
"boxmdary  conditions  which  assert  tnat  at  a  surface  of  separation  of  two  different 

homogeneous  media,  the  components  of  S  and  H  tangential  to  the  surface  shall  vary 

t 
continuously  from  one  side  of  the  surface  to  the  other  .   Since  our  boundary  sur- 
faces are  horizontal  planes  and  since  from  (8)  and  (9),  the  only  non-vanishing 
horizontal  components  of  3  and  H  are  S,  and  H/^,  it  is  to  these  components  that  i/e 

apply  the  conditions.   Let  the  separating  surface  lie  betv/een  the  j   and  (j  "*•  l) 

Xiql ^ — — 

^'Stratton,  op.  cit.  pp.  U85,  US6. 

'\ 

It  v/ill  be  noted  that  v/hereas  in  article  2  v/e  employed  a  boundary  condition  re- 
quiring the  continuity  of  the  normal  component  of  vW-H,  v/e  now,  for  convenience,  use 
an  eouivalent  condition  requiring  the  continuity  of  the  tangential  component  of  H. 
The  equivalence  of  these  two  conditions  is  proved  in  Stratton,  op. cit ., pp. 6,  3U-37. 


",    1. 
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layers  and  let   the  subscripts    (j)   and(,j  ■•■  1)   denote  the  evaluation  of  quantities 
in  the  respective  layers.      Then  on  the  boundary  surface,   we   require  that 

^r(j)     "  \(j  +  1)      •  ^^5) 

But  from  (8)  and  (9)  ^te   see  that  these  imply 

d)  z  3r         ^  z  dr 


(16) 


(17) 


^!      ^  ^J  .   "?-  1       ^  ^j  -M 


iGJU.    c>r       iCJ/x.^        ar 


(18) 


Equations  (I7)  and  (18)  constitute  the  boundary  conditions  on  ^^l  at  any  surface 
separating  two  of  the  strata.   They  may,  hov/ever,  be  simplified.   In  the  first 
place,  both  (I7)  and  (18)  hold  at  all  points  of  the  boundary  surface,  which  is  to 
say  that  they  are  true  for  all  values  of  r.   Hence  we  may  integrate  them  v/ith  re- 
spect to  r,  and  since  from  (I3),  tiie  function  V|/ vanishes  at  infinity,  the  constant 
of  integration  in  (IS)  is  zero.   Furthermore,  from  (lU)  and  (I3)  together  we  have 
the  fact  that  -r^    vanishes  at  infinity  and  hence  so  does  -zr-!-  .   Thus,  the  con- 

stant  of  integration  in  (IJ)  is  also  zero.   Secondly,  we  set  the  magnetic  permeabil- 
ityyA  equal  to  unity  in  all  the  layers,  corresponding  to  the  physical  situation 
which  is  almost  invariably  encountered.  The  introduction  of  these  simplifications 
into  (17)  and  (18)  then  leads  to  the  following: 

c)  z  a  z        »  ^^^^ 

^J  "^j  "  ^j  ^  1  ^J  +  1  •  ^^  ^   0.1.2....  n)      (20) 

Equations  (I9)  and  (20)  represent  the  final  form  of  the  boundary  conditions 
v/hich  our  solution  function  rmst  satisfy  at  the  separation  surface  of  any  tx/o  of  the 
strata. 

By  v;ay  of  recapitulation,  the  mathematical  statement  of  the  propagation 
problem  of  Fig.  1  is  the  folloi/ing.   V/e  seek  a  function  ^V  which  satisfies' conditions 


A,B.C,  and  D  described  atove,  or  more  specifically,  v/hich  satisfies  equations  (lO), 
(12),  (13),  (lU),  (19),  and  (20),   If  we  can  find  such  a  function,  then  the  field 
vectors  E  and  H  derived  from  it  according  to  the  formulas  (S)  and  (9)  or  (5)  and 
(6)  constitute  the  solution  to  our  propagation  problem. 

It  uill  be  convenient  to  deal  v/ith  this  problem  by  dividing  it  into  two 
main  cases.   In  the  first  case,  v/e  consider  the  problem  of  Fig.  1  v/ith  the  dipole 
in  the  homogeneous  part  of  the  atmosphere  and  therefore  outside  of  the  m\iltilayer 
region.  This  case  is  of  importance  to  the  problem  of  airplane  to  ground  communica- 
tion.  In  the  second  case,  we  consider  the  problem  of  Fig.  1  with  the  dipole 
located  anyv/here  v;ithin  the  multilayer  region.  Section  II  will  be  devoted  to  the 
solution  of  the  first  of  these  cases. and  the  second  case  will  then  be  solved  in 
Section  III. 

Section  II 

Solution  of  the  Propagation  Problem  with  the  Dipole  Situated 
outside  the  Multilayer  Region. 

The  basic  procedure  associated  v/ith  the  method  we  shall  use  consists  of 
two  parts.  First  we  shall  solve  a  plane  wave  problem  v;hich  is  much  simpler  than 
that  of  the  propagation  from  a  dipole  source.  The  solution  of  the  problem  of 
propagation  of  dipole  radiation  is  then  "built  up"  out  of  the  solutions  of  the 
simpler  problem.   In  the  first  part  of  this  section  ue  shall  describe  and  solve 
the  plane  v/ave  problem.  The  superposition  principle  by  which  the  solution  of  the 
dipole  problem  is  built  up  is  then  explained,  and  the  solution  is  obtained  for  the 
case  in  which  the  dipole  is  situated  outside  the  multilayer  region, 

6,  Descri-ption  of  the  Plane  "Jave  Problem. 

Referring  to  Fig,  2,  0  is  the  origin  of  an  x.y, z-coordinate  system  with 
the  positive  y-direction  into  the  paper.   The  negative  z-region,  i.e.  the  region  of 
space  belov;  the  x;''-plane,  is  filled  with  a  homogeneous  medium  extending  indefinitely 
in  the  x-,y-,  and  negative  z-directions  and  having  electromagnetic  constants  £,  ^ 
<S      (the  magnetic  permeability  U,    is  taken  as  unity  throughout.)  The  positive 
z-region  consists  of  n  layers  or  slabs  extending  indefinitely  in  the  x-  and  y- 
directions  and  having  thicknesses  d, ,  d„,  ...  d  in  the  z-direction.   The  boundary 

surfaces  of  these  layers  are,  of  course,  planes  parallel  to  the  xy-plane,  and  the 


% 


0\'j 


15. 


i. 


n 


t; 


fr 


"^ 


li. 


.i^L_ 


vT 


A^ 


e.,^. 


Fig.    2 


J^' 


-th 


electromagnetic   constants  of   the  j         layer  are    g.,   C.    .     i^ljove  this   n-layer,    or 

multilayer,    system  is  another   semi-infinite  homogeneous  mediiom  rnth  constants 

£  1  .    'S'      .    T    .      Our  configuration  thus   comprises   a  multilayer  system  of  media 

n  ■•"  1         n  ■*"  1 

sandv/iched  "between  tv;o   semi-infiuite,   homogeneous   media. 

Let  us   now  suppose   that   from  the  negative   z-region  a  continuous   stream  of 
plane   electromagnetic   \raves   of  constiint   amplitude  and  polarized  so   that   the  H  vector 
is   always  parallel  to   the  xy-plane,    advances   toward  the   multilayer,    impinging  upon 
the  latter  at   the   surface   z  =   0.      These  plane  waves   all  have  parallel  wave  fronts, 
and  their  com::ion  normal  maJces   an  angle       Q     v;ith  the   z-direction  rs   shown   in  Fig   2, 


where  5  S'    represents   a  t:,"pical  \7ave   front    and 


iJ  tne   nor;.ial   to   it. 
..0 


'The  pro- 


jection of  xi  upon  the  xy-plane      is   also  presumed  to  make   an  angle  (J)     \;ith  the 
positive  X-axis,    although  this   angle   is   not    shovm  in  Fig.    2.        Q     and    (jJ     are  thus 
the   ordinary   spherical  polar   angles   of  the    direction  K,    these    Deing  su.fficient    to 
specify  the  direction  uniquely.      The   superscript  0  is  used  to   denote   these   angles 
in  the   region  z   f   Zq   .    (ilote  that    in  Fig.    2   z^  =   0.  ) 

The  Waves,    as   they  pass   through   the  iiw.ltilayer ,    will  undergo   reflection 
and  refraction   at   each  interface   so   that    the   disturliance   at   any  point    (x,y,z),   at   a 
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given  time  t,  will  consist  of  the  superposition  of  the  disturbances  due  to  all  the 
reflected  and  transmitted  waves  passing  through  (x,y,z)  at  the  time  t.   Since  the 
incident  waves  impinge  on  the  multilayer  in  a  coiitinuous  and  uninterrupted  manner 
with  constant  amplitude,  and  since  at  a  fixed  point  in  space  the  incident  disturTj- 
ance  is  harmonic  in  time,  ue  may  assert  that  the  disturtance  anywhere  v/ithin  or  out- 
side the  multilayer  \/ill  ultimately  settle  down  into  a  steady  state  wherein  the 
disturbance  at  any  point  is  harmonic  in  time  and  has  an  amplitude  that  is  a  function 
of  position.  The  truth  of  this  assertion  is  intuitively  evident  on  physical  grounds. 

We  may  now  nsk  the  question:  given  the  configuration  shown  in  Fig.  2^ 
\/hat  is  the  amplitude  of  the  ultimate  disturbance  as  a  function  of  position,  i.e., 
what  are  the  forms  of  the  vector  functions  E(x,y,z)  and  H(x,y,z)  in  the  expressions 

Ee      and  He      vi/hich  constitute  the  ultimate  disturbance  under  consideration? 
This  is  the  plane  \rave  problem  v/hich  \ie   set  ourselves,  and  its  importance  resides  in 
the  fact  that  from  its  solution  that  of  the  dipole  propagation  problem  can  be  con- 
structed. 

This  simpler  problem  has  been  completely  solved  by  R.  K.  Luneberg,  and 
his  solution  is  described  in  another  report  devoted  entirely  to  this  subject. 
Research  Report  Ko.  I72-3.   Indeed,  Limeberg  deals  with  a  somewhat  more  general 
problem  than  the  one  of  Fig.  2   in  that  he  does  not  restrict  the  incident  plane 
i/aves  to  those  polarized  './ith  their  H  vectors  parallel  to  the  xy-planc ,  but  admits 
incident  plane  waves  of  arbitrary  polarization.  Also  *-  although  this  point  is  not 
of  great  significance  -  he  allows  the  magnetic  permeability  to  differ  from  layer  to 
layer,  \;hile  we  have  assumed  it  to  be  unity  throughout. 

Because  of  the  arbitrary  polarization  of  the  incident  plane  v/aves  in 
Luneberg's  treatment,  it  is  not  possible  for  him  to  describe  .the  solution  of  his 
problem  in  terms  of  a  single  scalar  potential  function.   It  is  for  this  reason  that 
we  have  restricted  the  problem  of  Fig.  2  to  the  case  v/here  the  H  vector  of  the  inci- 
dent wave  is  parallel  to  the  xy-plane.  Under  such  a  restriction  our  incident  field 
is  transverse  magnetic  lath  respect  to  the  z-direction,  and  hence,  by  the  reasoning 
already  outlined  in  the  first  section,  the  solution  field  of  the  plane  v/ave  problem 
is  also  transverse  magnetic*   But  \/e  knov/  that  this  transverse  magnetic  property 
implies  that  there  exists  a  scalar  potential  function  -  the  Hertz  potential  -  from 
v/hich  the  solution  vectors  E  and  H  of  our  plane  \/ave  problem  can  be  derived  by  means 
of  formulas  (5)  and  (6)  or  (2)  and  (3).   Thus  our  plane  wave  jjroblem  has  been  re- 
duced to  the  finding  of  a  single  unknovm  scalar  function,  exactly  as  was  the  general 
propagation  problem  of  Fig.  1. 
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7.  iHathematical  Expression  for  Plane  ^/aves. 

Before  going  further  \/ith  the  plane  v;ave  problem,  let  us  write  a  mati^e- 
matical  expression  which  will  represent  the  incident  plane  wave  stream,  or,  for 
that  matter,  any  plane  v;ave  stream.   If  (^  and  ^  denote  vectors  \;hich  are  condtaxit 
in  both  space  and  in  time,  then  the  vector  form  of  the  incident  v/ave  system  v.'ill  be 

_*   -^^t     i5  lkj^(x  sin  Q  cos  <?*  +  y  sin  0  sin  CP     +  z  cos  G  )  -i  oJt 

E  e"'"'     =  <!•  e  (21) 

.    .  ^,     —  ik  (x  sin  G  cos  (0     +  y  sin  G  sin  (fi     +  z   cos  0  )  -iCJt 
H*  g-iC^t  ^  ^   0 

e 

From  the  form  of  the  exponential  we  see  that  (21)  represents  a  disturb- 
ance whose  constant  phase  surfaces,  or  wavefronts,  are  planes  the  common  normal  to 
which  has  a  direction  specified  by  the  polar  angles  0  and  (p    .     Moreover,  by 

suitably  relating  the  amplitude  vectors  ^  and  ^,   it  is  easy  to  make  E  e~ 

— *  -i<Ot 
and  H  e      satisfy  haxwell's  equations.   Thus  equations  (21)  are  of  the  correct 

form  for  representing  the  incident  plane  \/aves.  V/e  desire,  however,  to  represent 

the  incident  wave  system  by  means  of  a  single  scalar  function,  rather  than  by 

several  vector  functions.   If  \;e  denote  such  a  scalar  function  by  g*,  then  evidently 

this  function  must  have  essentially  the  same  form  as  E  and  H  ,  i.e., 

-icJt    0  ^^0^^  ^^^  ®  ^°^  ^     "^  ^  ®^^  ®  ^^"^  9?  +  z  cos  0  )-icjt 
g*  e~     =  A  e  .        (22) 

where  A  is  a  scalar  amplitude,  constant  in  both  space  and  time. 

Application  of  forraulas  (5)  and  (6)  to  equation  (22)  leads  directly  to 
expressions  for  S  and  H  that  are  identical  in  form  with  those  of  (21 ),  and  which 
automatically  satisfy  i'iax\;ell's  equations.   From  (22)  we  see  that  the  incident  dis- 
turbance, even  v/hen  expressed  in  x.erms  of  the  Hertz  potential  still  consists  of  a 
stream  of  plane  waves.   Thus,  our  reformulation  of  the  plane  wave  problem  in  terms 
of  a  scalar  potential  function  has  not  altered  its  plane  wave  character;   this 
character  ivill  become  tiore  evident  in  \/hat  follows. 

8.  Solution  of  the  Plane  ',Jave  Problem. 


Let  us  nov/  proceed  directly  to  the  solution  of  the  plane  wave  problem. 
Instead  of  stating  this  problem  formally,  hov/ever,  as  \;as  done  in  the  case  of  the 
propagation  problem  in  Section  I,  we  shall  anticipate  natters  by  assuming  the  form 
of  the  solution  function  and  shov/ing  that  this  function  can  be  made  to  satisfy  all 
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the  conditions  uhich  the  protlen  inposes  upon  it.     For  convenience,   moreover,   \;e 
take   the   incident  v/aves   to  be    of  unit   amplitude,    i.e.,    in  equation   (22),   A     =   1. 
With  this   assumption  let  us   denote  the  solution  of  the  plane  v;ave  problem     by 
g  e  ,   \;nere  g  may   oe  v/ntten 

g  =  g(x,y,z;      0°,    ^°)  (23) 

Equation   (23)    simply  expresses   the    fact   that    the   amplitude   of  the  ulti- 
mate  disturbance   at   any  point   uithin  or  o\itside  of  the  j.Tultilayer   depends  upon  the 
coordinates    (x,y,z)   of  the  point   in  question  and  also  upon  the   direction  angles 
e   ,    (y     of  the   normal   to  the  plane  v/ave  fronts   of  the   incident  \7aves.      The  choice 
of  a  particular  fonn   01    the  function  g  \;ill  be  based  on  the  physically   reasonable 
expectation  that,   uhen  a  plane  v/ave   impinges  upon  a  multilayer  system  such  as  that 
of  ?ig.    2,    the  ultimate   disturbance  in  any   layer  v/ill  be  made  up  of  tv/o  plane  \mves; 
a  forward  moving   or  transmitted  wave  and  a  back\/ard  moving  or  reflected  wave.     We 
should  also  expect  that,    if  in  the  j       layer  the  normal   to   the   forward  moving  plane 
wavefront  makes   an  angle  G     with  the  positive   z-direction,    then  in  accordance  with 
the  ordinary   reflection  la\/,    the  angle  betv;een  the  normal   to  the  refl'ected  \/avefront 
and  the  positive   z-direction  will  be  tt  -  G*^.      Finally,   we   should   expect    that,   v/hile 
the  angles  9  of   the  various  wavefront   normals  v;ill   suffer  a  change  throu^  refrac-  ■ 
tion  as  we  cross   the  boundary   surface  separating  tv/o  layers   (i.e.,   0     j-  0  ), 

the  angles  (^'of    the  same  normals  \/ill  all  be  equal  to    ^   .      The   last  assxunption 
means  that   the  various   reflected  and  refracted  wavefront  normals   or  rays,   as  we  may 
call  them,   all  lie  in  the  plane   of  incidence   of   the  original   incident  wave. 

All  of  the  foregoing  assumptions  will  be  justified  a  posteriori.      They 
v/ill  now  be  given  mathematical   expression  through  the  assertion  that  they  are  equi- 
valent  to   the   reouirement   that    the  ultimate   disturbance   have   the   form 

.,^  .     ik  (x  sin  0*^   cos  <^      +  y   sin  ©"^   sin   <^     +  z  cos  Q"^ )   -i6Jt 

g  e~  =  A*^   e     "^ 

.      ik  (x  sin  0"^   cos  ^     +  y  sin  0"^   sin  Up  -  z  cos  ©"^ )   -i6Jt        (2U) 
+     B"^    e     "^ 

(j   =  0,1,2,    ...    n  +  1) 

It  will  easily  be  seen  by  inspection  that  (2U)  embodies  all  of  the 
assumptions  described  above.   The  first  term  on  the  right  represents  the  forward 
moving  v/ave  in  the  j   layer,  and  the  second  term  the  corresponding  reflected  v/ave. 
In  the  semi-infinite  negative  z-region  (j  =  0)  the  first  term  on  the  right  represents 
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0 
the  original  incident  wave  for  \/hich  by  assumption  A  =  1.   In  the  other  semi-infin- 
ite medixim  (j  =  n  +  l)  there  will  he  no  reflected  wave  since  there  is  no  upper 

boundary  surface,  and  hence  3     =0.  The  remaining  2n+2  unkno\;n  amplitude  coeffi- 
cients A  and  5  are  functions  of  9  and  (p   ,  and  \;e  shall  nov/  show  ho\/  they  are 
determined  by  the  2n+2  boimdary  conditions  at  the  various  separation  surfaces  of 

-io* 
our  multilayer  systen;  so  that  g  e     \i\\\   be  the  solution  to  the  problem  of  Fig.  2. 


Let  us  suppress  the  tine  factor  in  (2U)  and  v;rite 

ik.(x  sin  Q-^cos  V^°+  y  sin  0*^  sin^'^  +  z  cos  9"^) 

LK.^x  sm  wcos  '       

B-^  e   ' 


g(x.y.z:  e°.^°)  =  A-^  e  ^J 


.   ik.(x  sin  Q-^cos '/>  +y  sin  9*^  sin(y?  -  z  cos  9*^)     (25) 
+  B"^  e   "^ 


We  shall  determine  the  coefficients  A  and  3*^  so  that  the  function  g 


satisfies  our  problem.  V/e  first  observe  the  following:  Denoting  by  g*^  the  values 
whic 


which  g  takes  on  in  the  j   layer  there  results  by  direct  substitution  in  equation 


V^  gJ  +  k^  gj  =  0.  (26) 

That   is,    regardless  of  the  coefficients  A     and  3*^,    g  satisfies   the    v/ave  equation 
(26)   as   required,    everyv/here   in  the   j        layer, 

For  the   determination  of  A     and  B     we  make   use   of  the  boundary  conditions 
at   the  various   separation  surfaces   of  the  multilayer.      These  boundary  conditions   are, 
as  stated  previously,    the   continuity  across  the  surfaces  of  the  horizontal   components 
of  the   electric  and  magnetic  fields.      3y  fomulas   (5)   and  (6)   the   horizontal 
components    in  question  are: 

^2_  ^2 


X  az  dx        •  ^y  dz  by        ' 

(27) 
,  2        ^  2 

X        icj     dy  '  y  ito     ax  • 

Now  from  the  assumed  plane  wave  nature  of  the  phenomenon,  it  is  evident 
that  the  dependence  of  the  amplitude  of  the  ultimate  disturbance  upon  position  in 
any  plane  parallel  to  the  plane  of  incidence  will  be  independent  of  the  orientation 
of  the  plane  of  incidence  with  respect  to  some  fixed  coordinate  system,  i.e.,  will 

be  independent  of  ^  .   Thus  the  unknovm  coefficients  A"^  and  3*^  are  functions  of 
9  alone,  rather  thaii  of  both  9  and  ^      as  was  stated  above.   For  simplicity  in 
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dealing  v.'ith  our  boundary  conditions,  therefore,  v/e  sliall  choose  Cp  =  0,  That  is, 
the  plane  of  incidence  is  the  xz-plane.  Since,  in  this  case  g  no  longer  depends  on 
y,  equations  (27)  "become 

?.  2 

\  ~      .^2  3x  \  ito   37  ^^^^ 

Application  of  the  "boundary  conditions  at  the  surface  separating  the  j    and 

st 
(j  +  1)   layers  nov/  leads  to 


^'g-^    _      ^V*' 


3z  (3x        J  z  3x 


(29) 


•^    dx         j  +  1      ax      •  ^^^'' 

and  integration  of  (29)  and  (30)  uith  respect  to  x  results  finally  in 


a 


r")  t::"        _   ag" 


a  z        3  z 


(31) 


^^^  =  ^^  ,  g^  *  ^  •  (32) 

(j  =  0,1,2  ...  n) 

Equations  ( 3I )  and  (32)  ex]Dress  the  fact  that  the  Hertz  potential  g, 

the 
which  constitutes  the  solution  of/plane  v/ave  pro"blem,  satisfies  the  same  "boundary 

conditions  at  the  various  horizontal  interfaces  as  does  the  Hertz  potential W  v/hich 

represents  the  solution  to  the  propagation  pro'olem  of  Fig.  1,  the  conditions  on  W 

"being  given  by  (19)  and  (20)  .   It  is  to  "be  noted,  incidentally,  that  (3I)  and  (32) 

hold  for  all  values  of  (p   ,    the  latter  having  "been  set  equal  to  zero  only  for  con- 
venience in  olotaining  our  result, 

'  In  integrating  (29)  and  (3O)  to  give  (3I)  and  (32),  \/e  have  set  the  constants  of 
integration  equal  to  zero.  This  procedure  may  be  justified  either  by  means  of  a 
physical  argument,  into  \fhich  \;e  slaall  not  .enter,  or  simply  by  observing  that  either 
(29)  and  (30),  or  (3I)  and  (32),  lead  to  the  saroe  equations  (33)  and  (3U)  if  account 
be  taken  of  condition  (35) • 
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Applying  conditions  (3I)  and  (52)  to  the  expression  (25)  for  g  at  each 
of  the  separation  surfaces,  v;e  have 

.   ik.(x  sin  9*^  cos^^  +  y  sin  9"^  sin  ^  +  z  .  cos  O*^) 

A''  ik.  cos  9"^  e  "^  »> 

J 

.  ik  (x  sin  9"^  cos  (/V  +  y  sin  9*^  sinCt*^  -  z  .  cos  O*^) 

-  3"^  ik.  cos  9=^  e  -^  J 

J 

.^i  .  ,   ik.   (x  sin  9"^"*"^  cos  (^  +  y  sin  9'^'^'^   sin'/^°  +  z.  cos  O'^'^-'-) 

=  A"^   ik.^,  cos  9*^   e  "^  J 

j+1 

..-,  ..-,   ik.^,  (x  sin  9*^   cos  (^'  +  y  sin  9^^  sinCf'  -  z.cos  O*^"*"  ^ 

-  B*^   ik.^^  cos  9"^   e  ^  J        (33) 

.      ik.(x  sin  0''  cos  1^   +  y  sin  0*^  sin  ^  +  z.  cos  0"^) 
A-^  k":  e  -^  ^ 

.      ik.(x  sin  9*^  cos  (//'^  +  y  sin  9"^  sin  <p^  ~  z  .   cos  9*^) 
+  B''  k7  e  -^  "^ 

•+1  2     ik.   (x  sin  9'^"^-'-  cos  ^4-'^  +   y  sin  9'^"*'^  sin  (p^  +  z  .   cos  e^'*''^) 
=  A*^   k._^,   e  "^  ^ 

ik.   (x  sin  9'^^-'-  cos  CP^  +  y   sin  9'^*-'-  sin  (f  -  z .   cos  e^"^""-) 
+  B-J   k^^^  e  -^  J  (3U) 

(j  =  0,1,2,...    n) 

In  these  expressions   z.   represents  the  height   atove  the  xy-plane  of  the   surface 

th         '^  s  t 

separating  the   j        and   (j+l)        la^'ers;    thus   z^  =  0. 

Equations    (33)   and    (3U)   constitute   Lhe   2n+2   relations   from  v;hich  the 
2n+2  unkno\-/ns  A      and  3      are   determined.      These   relations   may  "be   simplified,   hovrever, 
in  view  of  the  fact   that    (33;   and  (3U)  must  hold  everywhere  on  the   interface 
z   =   z.,    i.e.,    they  ciust   hold   for  all  values   of  x  and  y  vrhen   z  =   z..      Hov/,    tv/o  such 
functions   of  x   and  y  as   those    on  the   left   and  right   sides    of   (33)   and   (314-)  can  "be 
equal   for  all  values   of  x  and  y   if,    and  only   if,    the   coefficients   of  x   and  y   in  the 
exponientials   on   one   side   are    equal   to   the   corresponding   coefficients   in  the   expon- 
entials  on  the   other   side,      ihus   v/e  :mist   have   ik .   sin  9*^   cos   (P    =ik     ,sin  9"'"^  cos  <y-' 

J  J+1  ^ 

and  ik.    sin  9      simi^'     =    ik  ._^,    sin  9*^        sin  ^/  ,    "booh   of  v/hich  reduce   to 

k.    sin  9''   =  k.        sin  9'^'^-'"      .  (35) 
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Equation  (35)  is,  of  course,  Snell's  refraction  law,  which  ue  now  see 
must  be  true  for  plane  waves  of  Hertz  potential  as  v/ell  as  for  ordinary  electro- 
magnetic \/aves  if  the  function  g,  given  by  (25),  is  to  solve  the  problem  of  Fig.  2. 

Substitution  of  (35)  into  (33)  and  (3U)  reduces  these  equations  to  the 
final  simplified  form: 


.  ik.z.cos  0*^   .       .  -ik.z.cos  9"^  ik .   z  .cos  0*^ 

A-^k.  cos  O-^e  "^  "^     -  B-^k.cos  Q-^e   "^  '^      =  A"^   k .  ,cos  0*^   e  ^        "^ 
J  J  J+1 

I 

i+1 

-  B-^  >-.  cos  9-^  ^e   -^   -^  (36) 

.  ^  ik.z.  cos  9*^    .   ^  -ik.z.  cos  9"^     .  -,  ^   ik .  ,  z .  cos  9*^ 
J  J  J"^l 

i+1 

-ik  .^,  z  .  cos  0*^ 


-   bJ"H:2  ^  e   "J"l  J 
J^l 


(37) 


(J  =  0,r,2,...n) 

Using  (36)  and  (37)  and  remembering  that  A  =  1,  B    =  0,  \/e  can  solve  for  all 

the  unknowns  A'^  and  B  ,   By  the  use  of  (35)  #  ii'^  which  the  values  of  k.  are  known  in 

terms  of  C .  ,    u . ,    and  C,  the  values  of  9"^can  be  expressed  in  terms  of  9  . 

9.  Summary  of  the  Plane  Wave  Problem. 

Let  us  s\iminp.rize  our  results  thus  far.   We  have  a  function  g(x,y,z;9,  6^/) 
defined  by  (25)  possessing  the  properties: 

(a)  V   g  *  k  g*^  =  0  everywhere  in  the  j   layer.  ' 

J 

(b)  If  A*^  and  B''  satisfy  (36)  and  (37),  then  g  satisfies  the  boundary  conditions 
(31)  and  (32), 

(c)  In  every  layer  g  represents  the  sum  of  a  refracted  plane  v/ave  and  a  reflected 
plane  wave.   In  particular,  in  the  semi-infinite  negative  z-region  (j=0),  g  com- 
prises the  sum  of  the  original  incident  plane  wave  of  unit  amplitude  and  a  reflected 
plane  wave. 

In  viev;  of  (a),  (b)  and  (c)  \ie   may  say  thaL  g  constitutes  the  solution 
to  our  plane  i/ave  problem  of  Fig.  2  expressed  in  terras  of  the  Hertz  potential.   V/c 
see  nov;  that  the  physical  assumptions  made  on  pages  I9  and  20  regarding  the  expected 
behavior  of  these  Hertz  potential  waves  were  all  justified.   They  obey  the  same  re- 
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flection  and  refraction  laus  as  do  light  waves,  at  least  for  configurations  such 
as  that  of  Fig.  2. 

Having  solved  the  plane  v/ave  problem,  \re   nov;  proceed  to  an  explanation 
of  the  superposition  principle  and  of  the  method  wherety  this  principle  is  utilized 
so  as  to  constmct  the  solution  of  the  dipole  problem  from  those  of  the  plane  v;ave 
problem, 

10.  The  Superposition  Principle. 

First  let  us  ooserve  that  "both  the  diffsrential  e'luation  (26)  and  the 
boundary  conditions  (3I)  '\nd  (32)  satisfied  by  g  are  linear.  From  this  fact  ve   may 
drav/  tv/o  conclusions: 

(a)  If  g,  and  gp  are  t\/o  different  functions  satisfying  the  differential  equation 
and  the  boundary  conditions,  then  g^  +  g  also  satisfies  the  differential  equation 
and  the  boundary  conditions, 

(b)  If  g  is  any  function  solving  the  plane  uave  problem  v/ith  incident  v;ave  of  unit 
amplitude,  then  A  g  solves  the  plane  wave  problem  uith  incident  \/ave  of  amplitude 
A°. 

It  is  clear  that  \;e  may  combine  (a)  and  (b)  to  arrive  at  the  follov/ing 
more  general  conclusion.  Referring  to  Fig,  2,  let  there  be  incident  upon  the  sur- 
face z  =  0  of  the  multilai'er  H  plane  waves  having  amplitudes  A^  ,  Ap,  . . .  A.^  and 
corresponding  propagation  directions  given  by  the  pairs  of  polar  angles 
(0^,  ^-j^),  (6p,  (PJ (9^^,(^^).  Then  we  may  say  that  the  form  of  the  ultimate 

disturbance  v;hen  these  K  [olane  waves  impinge  on   the  lovfer  surface  of  the  multilayer 
from  the  negative  z-region  is 


^=  aJ  g(x.y.z;  e°.(^°)  +  A°  g(x.y,z;  Q^,  (p^)   +  . . .  Aj°  g(x,y,z;  0°,  <f)^)   .        (38) 

That  is,  the  ultimate  disturbance  is  simply  the  algebraic  sum  of  the  separate  ulti- 
mate disturbances  due  to  each  of  the  incident  plane  waves  acting  separately;   in 
short,  superposition  is  valid. 

The  significance  of  this  result  is  apparent.   It  means  that  having  solved 
the  plane  v/ave  problem  in  v/hich  the  incident  field  is  a  single  plane  wave  of  unit 
amplitude,  v/e  can  immediately  write  down  the  solution  to  the  more  general  problem 
v/here  the  incident  field  is  any  field  which  can  be  expressed  as  the  sum  of  plane 
\/aves  of  arbitrary  amplitudes  and  directions  of  propagation.   It  is  necessary  to 
know  only  the  f^jnction  g(x,y,z;  0  ,  (f),    i.e.,  to  know  the  coefficients  A*^  and  B"^ 
of  equation  (25). 
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V/ith  the  validity  of  the  superposition  principle  established  for  the 
case  uhere  the  incident  field  is  representalile  as  the  sum  of  only  a  finite  num'ber 
of  plane  v/aves,  it  is  natural  to  ask  whether  the  principle  will  still  "be  valid 
for  the  case  of  a  more  general  type  of  incident  field.   In  particular,  \;e  shall  "be 
interested  in  incident  fields  in  uhich  the  finite  sum  of  plane  \;aves  is  replaced 
by  a  certain  kind  of  definite  integral.   Deferring  considerations  of  the  validity 
of  such  a  replacement  for  the  moment,  let  us  see  ho\/  this  transition  is  effected 
and  v;hat  its  consequences  will  be  should  supejrposition  still  prove  to  be  valid. 

ViTe  first  write  the  expression  for  the  incident  field  consisting  of  a 
finite  number  of  plane  waves,  i.e., 

1,1^   v —   p.  ikpCx  sin  9  cos  '   +  y  sin  ©  sini/''  +  z  cos  6  ) 

^     =  l__    k     e  s/s  s  s  s  ^^^^ 

s=l       ^ 
In  this  expression  the  summation  extends   over  all  the   discrete  values   of  ©   ,    C^ 
and  A     v;hich  correspond  to  the  integral  values   of  s  from  1   to  IT,      If,    instead  of 
permitting  these  variables   to   take  on  only  these  discrete  values,  we  allov/  6     and 

r       to  vary  independently  and  continuously  over  the  region  ©     <,  Q     ^  ©^    . 

a     -  b 

T^   -  </^  -  (P,     ,  and  consider  A  to  be  a  kno\m  function  of  ©  and  </*  ,  (which  will 

usually  be  analytic),  the  expression  analogous  to  (Ul)  is 


tii*    /     "^       o/^^   n   n  ik  (x  sin  ©°cos(^+  y  sin  ©°sin(i«^+  z  cos  e°)  ^ 

^      =/    d(p7   A(©°.  ^°)e  °  ^  ^  de°  (U2) 

ra       a 

The  incident  field  of  (1+2)  requires  the  sum  of  a  continuous  infinity  of  plane  v/aves. 
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Referring  once  ,.:ore  to  Fig.  2,  suppose  that  incident  upon  the  lower 
surface  of  the  multilayer  is  a  field  ^)   of  the  type  given  "by  (U2).   If  we  assume 
that  superposition  is  still  valid  in  this  case,  then  the  ultimate  disturbance  il)  , 
by  analogy  \;ith  (3S),  \/ould  be 


A(©°.  '/^^)     gCx.y.z;  9°,  ^°)  d0°         {kj>) 


i.e.,    (U3)   v.'ould  "be  the  solution  of  the  problem  of  Fig.   2  uith  incident  field  given 
by  (U2). 

One  further  generalization  v/ill  be  needed.      Consider  the    integral   repre- 
sentation  of  the  incident    field  ly     as   given  by   (U2).      In  this  expression  the 

variables  of  integration  9     and.  (J)     take  on  all  real  values  between  the  real  limits 

of  integration;      the  paths   of  integration  of  both  9     and  (p      in  (U2)   are  segments 

along  the  real  0     and  real  (^      axes.      Suppose,   however,    that  we  have  an  incident 
field  Y      permitting  an  integral   representation  of  the   saiie  form   as    (^-2)  but   re- 
quiring that    in  the  course   of  the   integration  the  variables  ©     and   (^     take  on  com- 
plex values   as  well  as   real   ones.      In  other  words,    suppose   that   in   (42)    the  paths 

of   integration,    instead  of  being  segnents  along,  the  real  9     and  iS)     axes,    are  con- 
tours  in  the   complex  9     and  U)     planes.     The  incident   radiation  would  then  be 


llJ*         I       c\  I         r\       cs     ^^n^^  ^^"^  ®  cos^  +  y  sin  9   sin^+  z  cos  9    )     ^ 
T      =    /    d^/   A(9".^'')e     "  ^  ^     ^  d©° 


iM) 


^^0  ^c^o 
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v;here  C  ^   and  C(77„  are  the  paths  of  integration  of  the  complex  variahles  0  and  P 
respectively,  Pig.  3»  ^^®  lengths  of  C  Qand  C^-^  being  either  finite  or  infinite. 


4>  /^^/^A/£r 


O    /^L.flf^£i 


Fig.  3 

Mq   nov;  ask  the  question:  If  in  the  configuration  of  Fig.  2  the  radiation 
incident  upon  the  multilayer  from  the  negative  z-region  is  given  "by  (UU),  what  is 

the  ultimate  disturbance,  ^"i      In  order  to  ansuer  this  question,  let  us  first  ob- 

* 
serve  that  the  incident  field  ^^  in  (UU)  is  no  longer  represented  as  the  super- 
position of  ordinary  plane  waves  since  the  exponential  form  in  the  integrand  of  (UU) 
now  contains  complex  angles  9  and  ^   instead  of  real  geometrical  ones,  and  these 
render  impossible  any  simple  physical  intferpretation  of  the  exponential  as  a  plane 
wave.  Hathem.-tically  speaking,  however,  the  situation  is  the  same  as  for  the  real 
case.   If,  in  Fig.  2,  we  go  back  to  the  original  problem,  where  a  single  plane  v/ave 
of  unit  amplitude  is  incident  upon  the  multilayer,  and  replace  this  problem  by  the 


.0 


0 


0  . 


,0, 


one  in  which  a  single  v/ave  of  the  type  e 


ik„(x  sin  0  cos*/   +  y  sin  0  sin  ij^^  +z  cos  G  ) 


is  incident  upon  the  multilayer,  0  and  '^-  now  being  complex  variables,  we  shall 
find  that  the  solution  of  this  problem  is  identical  in  form  with  the  solution  ob- 
tained when  0  and  ^  were  real.  The  ne\/  solution  v/ill  still  be  given  by 

gC^iy.z;  ©  t  ^'  )  of"  (25)  except  that  0  and  (p     have  the  appropriate  complex  values. 

It  i/ill  still  satisfy  the  same  differential  equation  and  the  same  boundary  conditions 

0       0 
as  it  did  when  0  and  (jp     were  real.   In  fact,  all  the  mathematical  properties  of 

g(x,y,z;  Q   ,<f-    ),   necessary  for  the  validity  of  the  superposition  principle 

ere  possessed  by  this  function  '.'hen  0   anci   C^  are-comple_x  .   .Even 
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Snell's  la\/,  equation  (35) »  carries  over  to  the  complex  case,  although  it  has  no 
apparent  geometrical  significance,  since  the  reasoning  used  to  establish  it  is  still 
valid. 

In  vievr  of  the  foregoing,  we  should  expect  that  if  the  generalization  of 
the  original  superposition  principle  to  incident  fields  representahle  "by  the  real 
integral  (U2)  is  justified,  then  its  generalization  to  the  case  v/here  the  incident 
field  is  expressed  "by  the  complex  integral  (UU)  is  also  justified.   If  this  expecta- 
tion is  correct,  ue  nay  surjnarize  our  conclusions  "by  saying  that  when  the  incident 
field  in  Fig.  2  is  given  by  {kh)   the  solution  field,  or  ultimate  disturbance  U-*  ,is 

It  is  clear  that  the  foregoing  discussion  is  contingent  on  the  validity 
of  the  superposition  principle  in  the  case  represented  by  equation  (44).  While  we 
have  shown  rigorously  that  this  principle  is  valid  v/hen  the  incident  v;ave  consists 
of  a  finite  number  of  real  or  complex  plane  waves,  v;e  have  not  as  yet  demonstrated 
mathematically  that  its  generalization  to  an  incident  v/ave  given  by  an  integral 
such  as  that  in  (UU)  is  justified.  This  demonstration  offers  no  particular  diffi- 
culty as  long  as  the  contours  C  ^  and  Ca-q  are  finite  in  length  and  pass  through 

y 

none  of  the  sing-alarities  of  the  integrand  in  (U5).  Under  these  conditions,  the 
possibility  of  the  integral  (U5)  or  its  derivatives  failing  to  exist  is  precluded. 
If,  however,  it  should  happen,  as  it  will  in  the  case  to  be  discussed  shortly, that 

one  of  the  contours,  say  C  ^  ,  must  go  to  infinity  in  some  direction  in  order  for 

0 

•* 
(l+U)  to  represent  the  incident  field  '-f-'  ,  then  (45)  becomes  an  infinite  integral 

and  the  question  of  its  convergence  and  that  of  its  derivatives  becomes  a  rather 
elaborate  mathematical  problem,  at  least  when  we  are  dealing  with  a  multilayer 
system  such  as  that  of  Fig.  2.  Unfortunately,  so  far  as  is  knovrn,  this  problem  has 
not  been  treated  in  the  case  of  the  multilayer.   Hence,  we  shall  confine  ourselves 
to  two  objectives.  (1)  V/e  shall  show  how  the  representations  (UU)  and  (U5),  if  their 
use  is  justified,  lead  to  the  solution  of  the  propagation  problem  when  the  source  is 
a  dipole  situated  outside  the  multilayer.  (2)  V7e  shall  make  more  explicit  the  par- 
ticular mathematical  difficulties  which  stand  in  the  way  of  justifying  this  general- 
ization of  the  superposition  principle. 


28. 
11,  Solution  of  the  Propagation  Problem  v:ith  the  Dipole  outside  the  Multilayer. 

For  convenience  \re    shall  locate  the  electric  dipole  in  the  homogeneous 
region  below  the  multilayer,  i.e.,  in  the  negative  z-region  of  Fig.  2.   It  is  assumed 
to  be  vertically  oriented  and,  vithout  loss  of  generality,  \;e  may  suppose  it  to  be 
placed  at  z  =  -h.   Ihe  radiation  from  this  dipole  impinges  upon  the  multilayer,  and 
our  problem  is  to  find  the  ultimate  disturbance  'y    anyv/here  v/ithin  or  v/ithout  the 
multilayer. 

In  order  to  solve  this  problem  Toy   the  use  of  the  generalized  superposi- 
tion principle,  we  must  first  \;rite  down  a  representation  of  the  type  (UU)  for  the 
incident  radiation  ^{-^    coning  from  the  dipole.  This  incident  or  primary  radiation 
is  simply  the  field  to  wliich  the  dipole  v/ould  give  rise  if  located  alone  in  a  medium 
of  constant.  k„  extending  indefinitely  in  all  directions.  In  the  first  section  of 

this  report  it  has  already  been  stated  that  this  field  is  given  by 

ik^R 

where  R  is  the  distance  from  the  obseirvation  point  to  the  point  (0,0, -h)  where  the 
dipole  is  located.  Equation  (U6)  therefore  gives  the  incident  radiation.   It  can 

^ik^R 
now  be  proved  that  the  integral  representation  (UU)  of  ^ is  given 

ty(2l)(22) 

2n       TT  -  1  00 


0  .'      q/  ik-     ik^h  cos  0  „   ik„(x  sin  0  cosCp  +  y   sin  G  sin^+  z  cos  9  ;    ^ 

(U7) 


®  -r^ =  /  dSi^l  --^^  e     "  sin  9"e      "  '  dfi 


for   z   >  -h  and  by 


Q     ik^h  cos  9  Q  ik^lx  sin  0  cosy?  +  y  sin  9  sin  ^'  +  z  cos  0  )   „ 

e  sin  9  e  39 

(Ug) 


for  z   <  -h. 


The  reason  for  requiring  one  representation  for  ^   v/hen  z  >  -h  and 


ik  R 

Einother  v/hen  z   <   -h  is   the   existence  of   the   singTilarity   in  0  at   R  =  0. 

"r 

j21)stratton,    op.    cit.   pp.    577,    578. 

^    '^'^Courant-Hilbert;    Methoden   der  ur-themat ischen  Physik,   Vol.11,   p.    I5U. 
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We  observe  that  in  both  representations  the  function  A  is  a  function  of 


0  alons,  and  is  the  same  in  both  (Uy)  and  (^-8),  narnely 

0  Ik  h  cos  G 

't/ith  regard  to  the  paths  of  integration,  ^'e  se 


(U9) 
e  that  C>^  is  the  same  in 


both  (U7)  and  (U8)  being  simpl/  the  segment  of  the  real/-'  -axis  from  0  to  2n.  Fig.  U 


sho\/s  C  p,  for  both  cases. 

e 


^  /='/,/y/\^^' 


/=/9r/v  /='o/^ 


Fig.  U 

We  are  no\i   in  a  ijosition  to  write  down  the  solution  of  our  problem. 
First  consider  the  region  z  >  -h.   According  to  (U5)  and  (U9)  the  solution  should  be 


S-' 


2n   ^  -  ioo  Q 

0  /  iko    ik^h  cos  0 
d(^  /   -^  e 


sin  9  g(x,y,z;  Q    ,  (^  )    dQ 


(50) 


lubstituting  for  g(x,y,z;  ©  ,^  )  from  (25)  v/here  A''  and  B"^  are  determined  "by   (36) 


-th 


^^fi  (37).  equation  (50)  gives  for  the  solution  in  the  j   slab  of  the  multilayer, 


'  d^I 
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TT 
,^_       —    -1  CX3 

2n     2 


„/ik_     ik  h  cos   9  n   -i     n     ik.(x  sin  ©''cos  (25  +  y  sin  ©'^sin(i^+z  cos  G"^)     ^ 


0        0 


(51) 

TT 


2n     2-^°° 


''^  Ak„     ik  h  cos  9  n    ^     n     ^'^A^  sin  Q'^cosW  +  y  sin  ©-^sin^i^-   z   cos  9-^) 

({//  -^  e     ^  sin  9^(9^)6     «^  ^  ^ 


0 


2n  -         dQ 

0  '0 
for  -  h  <  z  <  0  v/hen  j  =  0,  and  for  z  ,  <  z  <  z.  \/hen  j  =  l,2,....n. 

Hote  that  while  a  variable-  9*^ appears  in  the  integrand  of  (5I),  it  is 
in  reality  a  function  of  the  variaole  of  integration  9  throiigh  the  Snell's  law 

relation  k  sin  9"^  =  k„  sin  9  .   If,  however,  ue   solve  this  equation  for  9*^  to 

.  k 

o"btain  9*^  =  f(9  )  =  arc  sin  C^-  sin  9  ),  we  observe  that  f(9  )  is  not  single-valued; 

J 
i.e.,  the  Snell's  lav:  relation  does  not  uniquely  determine- 0*' as  a  function  of  the 
veiriable  of  integration.  To  remove  this  ambiguity,  the  following  rule  is  to  be 
used:  v/henever  the  real  part  of  the  variable  of  integration  lies  between  0  and  2,^ 

the  real  part  of  9*^  also  lies  between  0  and  -  ;  whenever  the  real  part  of  the 

variable  of  integration  lies  betv;een  ~  and  tt,  the  real  part  of  9*^  lies  between^ 

and  TT.   In  the  present  section  of  this  report,  the  real  part  of  9  lies  v/ithin  the 
range  0  to  ^  ,  so  that  only  the  first  part  of  the  rule  is  relevant.   In  Section 

III,  however,  ue  shall  see  that  the  real  part  of  the  variable  of  integration  may 
lie  in  the  range  ^  to  rr,  so  that  the  second  part  of  the  rule  must  be  used.  At  no 
time,  however,  will  the  real  part  of  the  variable  of  integration  lie  outside  the 
interval  0  to  tt.   With  this  rule,  we  are  assured  that  the  integrand  of  (5I)  is  a 
single  valued  function  of  9  and  6p  . 

We  have  therefore  found  the  ultimate  disturbance  any\/here  in  the  region 
z  >  -h  of  ?ig.  2  due  to  a  dipole  located  at  (0,0, -h),  assuming  of  co\arse  that  the 
superposition  principle  is  valid.   The  disturbance  is  given  by  (5I)  \/here  A*'  and  B'^ 
are  the  solutions  of  the  system  of  linear  equations  (36)  and  (37 ) .In  particular, 
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let  us  \/rite  the   expression  for  the  disturbance   in  the  region  "between  the  dipole  and 
the  multila^ver,    i.e.,   0   >   z   >   -h.      In  this   region  j   =   0  and  A*^   =  A     =1,    Then   (5I) 
becomes 


'ik^      ^^'-rP-  ^°^  Q  0   ^^0^'''  ^^^  ®   cos(^  +  y   sin  9   sin^  +  z   cos   6   )        ^ 

■— r  e  sin  6  c  d0 

2n 


(52) 


TT 

—   -     1  00 

2u  2 


ik^     ik_h  cos  9  n  />     o  i'^^(x  sin  9  cosW    +  y  sin  9  Rinip-  z  cos  «)    ^ 

+/d^/-^e      °  sinGV(9°)e°  ^  ^  d9°. 

(0   >   z    >  -h) 

If  we  examine  (52)  we  see  that,  according  to  (Uj),  the  first  integral 
in  (52)  is  merely  ikgR      Ecnce  we  may  write  (52): 

6  • 


R 


ik^R 
4^0  =   ^ 


TT 


.0 


„  ,  ik       ik  h  COG  9"         a  a     r     ik^fccsin  0  cos®+y  sin  9  sinftT-z  cos  9   )     _ 
^/d^7^e     °  sin  9^(9*^)  e     ^  ^  ^  d©° 

(53) 


The  interpretation  of  (53)  is  evident.   It  means  that  in  the  region  be- 
tween the  dipole  and  the  multilayer,  the  disturbance  consists  of  the  primary  dipole 
radiation  v/hich  v;ould  exist  if  there  were  no  multilayer  plus  a  "reflected  field"  due 
to  the  multilayer  and  which  is  represented  in  (53)  by  the  integral  term.  Physically, 

of  course,  this  division  of  the  field  into  a  primary  or  incident  field  and  a  re- 

ik  R 
fleeted  field  is  what  we  should  1-iave  exoected.   Moreover,  the  term     '0 

e 


R 
represents   the  primary  radiation  not   only   for   z    >  -h  but   also   for   z   <  -h,    and  since 
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we  should  expect  that  the  reflected  tenn  in  (53)  • 'together  with  its  derivatives,  is 

continuous  across  the  z  =  -  h  plane  / ^.•„      e  +  „   t  ^u   j-   >i    •   -,      ■  ^ 

[tne  - — 5 term  gives  the  dipole  singularity 

and  there  are  to  be  no  other  singularities,  there  being  no  other  sources)  ve   should 
conclude  t/iat  (53)  represents  the  disturbance  everyv;here  in  the  negative  z-region, 
and  is  not  restricted  to  the  region  betv/een  the  dipole  and  the  multilayer.  If  this 
expectation  is  correct,  then  (51)  sj^d.  (53)  together  constitute  the  complete  solution 
of  the  propagation  problem  of  the  dipole  situated  outside  the  multilayer.   The  ques- 
tion of  the  correctness  of  this  expectation  v/ill  be  found  to  be  a  part  of  the 
general  question  of  tne  validity  of  the  superposition  principle  represented  by  (UU) 
and  (U5),  for  the  case  i/here  one  of  the  contours  of  integration  goes  to  infinity, 
¥e  shall  now  discuss  this  general  question  briefly, 

12.  Validity  of  the  G-eneralized  Superposition  Principle. 

At  the  outset  let  us  say  what  v/e  mean  by  the  validity  of  the  superposition 
principle  represented  >y  expressions  (1+U)  and  (U5).  We  mean  simply  that  this  prin- 
ciple is  valid  if  any  solution  obtained  by  means  of  it  satisfies  conditions  A,  3,  C, 
and  D  of  Section  I  v/hen  the  antenna  is  a  single  vertical  electric  dipole.   Let  us 
examine  the  solution  Vl)  given  by  (5I)  and  (53)  from  this  point  of  view. 

th 
Condition  A  -  This  requires  that  in  the  j   layer  of  Fig.  2,  W    should 

2        2 

satisfy  \/     >^   +  k.  {O .  =   0.   From  (5I)  we  see  that  ^;  will  satisfy  this  equation  if 

it  is  permissible  to  perform  the  operation  \7  under  the  integral  sign,  sinoe  the 
integrand  itself  automatically  satisfies  the  equation  in  question.   Since  the  inte- 
grals are  improper,  the  contour  going  to  infinity,  this  interchange  of  the  operations 
of  differentiation  and  integration  CMinot  be  made  without  considering  the  convergence 
of  the  integrals ,   iJow,  in  the  first  place  \;e  do  not  know  tliat  the  integrals  in  (5I) 
are  themselves  convergent.   In  the  second  place,  even  if  these  integrals  are  con- 
vergent ,  v/e  do  not  kno\/  that  the  integrals  v/hich  result  from  performing  the  \/ 
operation  under  the  integral  sign  are  convergent.  In  order  to  resolve  these  converg- 
ence difficulties  it  is  necessary  to  have  a  much  more  precise  knowledge  of  the 
properties  of  the  fvmctions  A'^(0  )  and  3*^(9  )  than  that  they  merely  satisfy  equations 
(36)  and  (37)1  since  it  is  obviously  on  the  analytic  behavior  of  A''  and  B*^  that  the 
convergence  questions  depend.  Unfortunately,  we  do  not  as  yet  possess  the  knov;ledge 
of  this  behavior  and  therefore  cannot  deal  rigorously  \/ith  the  convergence  problem. 
This  difficulty  './ill  be  seen  to  arise  in  connection  with  the  remaining  conditions 
which  vy  is  supposed  to  satisfy. 
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Condition  B  -  This  condition. xeguires  that    in  the  inmediate  vicinity  of 

the  dipole   the   function  U/ "behave  like     e ,     According  to   (53).   this  condition 

R 
v/ould  evidently  mean  that  the   integral   term  should  remain  finite  and  continuous  at 

R  =  0,   a  condition  ue  have  already''  anticipated  in  connection  with  the   assertion  that. 

(53)  represents   the  solution  field  everywhere  in  the  negative  z~region.      Once  again. 


however,   we  see  that  the  convergence  of  the  integral  can  te   investigated  only 

through  a  knowledge  of  the  properties  of  the  function  3  (9    ). 

lim    'X' 
Condition  C  -  This   states  that  U.' must  satisfy     .^     ^    1     =  0  and 

'  K  ->  00 


lira 

R->  00  ^ 


^^--H' 


=  0,  Considerations  similar  to  those  explained  in  A  and  B 


evidently  apply  here  also. 

Conditions  D  -  These  are  the  "boundary  conditions  (I9)  and  (20)  vAiich  \U 
must  satisfy,  vre  have  already  seen  from  equations  (jl)  and  (32)  that  the  plane  waves 
satisfy  conditions  (19)  and  (20).  The  integrands  of  (51)  therefore  satisfy  these 
conditions,  and  hence  also  the  integrals,  subject  to  the  convergence  considerations 
already  discussed, 

13.  Conclusion  to  Section  II. 

In  this  section  we  have  used  the  method  of  superposition  of  plane  v;aves 
to  solve  the  propagation  problem  of  the  dLpole  located  outside  a  multilayer  region. 
The  solution  is  given  "by  expressions  (5I)  and  (53),  and  the  solution  vectors  f  and  H 
are  derived  from  expressions  (5I)  and  (53)  througli  the  use  of  the  formulas  (5)  and 
(6)  or  (g)  and  (9). 
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Section  III 

Solution  of  the  Propagation  Problem  v/ith  the  Dipole  Situated 
withija  the  Multilayer  Region 

Ik,    Statement  of  Pro"blem. 

In  this  section  v.'e  shall  be  concerned  with  the  solution  of  the  second 
case  of  the  propagation  prohlem  described  in  the  last  paragraph  of  Section  I.  Re- 
ferring to  Fig.  5 1  wi'-ich  is  merely  Fig.  1  with  added  details,  ve   have  0  as  the  origin 
of  a  rectangular  coordinate  system  \.'ith  the  positive  senco  of  the  y-axis  into  the 
paper. 


G  G  represents  the  surface  of  the  earth,  v/hich  is  also  the  xy  plane.  The 
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th 


vertical  electric  dipole  of  angular  frequency  fjJ  is  located  in  the  s   atmospheric 
layer  at  a  height  h  above  the  earth  and  is  situated  on  the  z-axis.   Our  problem  is 
to  find  the  ultimate  disturbanceV|^  anyv/here  in  any  of  the  layers. 

The  difference  between  the  present  problem  and  that  of  Section  II  is  that 
in  the  present  problem  the  dipole  is  located  v/ithin  the  multilayer  region  instead  of 
being  either  below  or  above  it.   V/hile  we  shall  find  that  the  general  method  developed 
in  Section  II  -  resolution  into  plane  v/aves  and  superposition  -  is  applicable  to  the 
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present  case  also^  the  new  location  of  the  dipole  requires  an  additional  procedure 
in  order  to  obtain  the  solution. 

15.  pecomposition  into  Subsidiary  Pro"blems« 

The  procedure  v;e  shall  employ  is  the  iollo\/ing.   V/e  brealc  up  our  propa- 
gation problem  into  tv/o  subsidiary  problems.   In  the  first  of  these,  only  the  upper 
half  of  the  dipole  is  considered  as  radiating;  i.e.,  referring  to  Fig.  5i  we  suppress 
the  primary  radiation  in  the  region  z  ^  ^  _  z  <  h,  considering  the  configuration  of 

Fig.  5  to  be  unaltered  other\;ise.   V/e  then  solve  the  problem  of  Fig.  5  under  these 
conditions.   In  the  second  subsidiary  problem  \;e  sujj^jress  the  primary  dipole  radia- 
tion in  the  region  h  <  z  _  z   and  assume  that  only  the  lower  half  of  the  dipole 
radiates.  The  sura  ci  the  solutions  of  these  two  subsidiary  problems  evidently  co)>- 
stitutes  the  solution  of  tne  problem  of  Fig.  5* 

Bach  01  these  subsidiary  problems  will  no\/  be  solved,  by  using  the 
method  of  Section  II  with  a  suitable  modification, 

16.  Solution  of  the  First  Subsidiary  Problem. 

In  order  to  solve  the  first  subsidiary  problem  through  the  use  of  our 
superposition  principle,  we  mast  first  find  the  form  of  the  plane  wave  function  g 
in  any  of  the  layers  in  Fig.  5«  ^^   begin  with  the  s   layer,  i.e.,  the  one  con- 
taining the  dipole. 

The  s   layer. 

In  that  part  of  the  s   layer  between  the  dipole  and  the  upper  boundary 
surface  of  the  layer,  i.e.,  h  i  z  £  z  ,  the  plane  wave  function  g  evidently  con- 
sists of  three  plane  waves: 

(a)  A  plane  \/ave  of  unit  amplitude  moving  in  the  same  sense  as  the  primary   . 
radiation.   This  will  be  identified  as  a  plane  wave  component  of  the  primary  radia- 
tion  v/hen,  in  the  final  integration,  v;e  multiply  g  by  the  function  A  of  equation 
(Ug)  of  Section  II. 

/  \  s  T    .     .    , 

(b)  A  plane  \:ave  of  unknovm  amplitude  B-.  ,  novmg  against  the  primary  radia- 
tion  and  representing  the  effect  of  reflection  from  all  the  layers  above  the  s 

^  •  'We  affix  the  subscript  1  to  the  amplitudes  associated  v/ith  the  first  subsidiary 
problem;  the  solution  to  the  latter  \/ill  be  called  \|>  .   Similarly,  we  shall  use  the 
subscript  2  v;hen  we  deal  v/ith  the  second  subsidiary  problem. 
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(c)  Another  plane  wave  inovir^  with  the  primary  radiation  and  of  unknown 

amplitude  A,  ,  representing  the  effect  of  reflection  from  all  the  layers  belov;  the 

th 
s  . 

The  function  g     can  thus  "be  written 


ik  (x  sin  ^     cos  W     "*"  y  sin  G     sinto    +  z  coc  0    ) 
g     =   (1  +  A^)   e  + 

ik  (x  sin  0     cos  ^    +  y  sin  9     sin(o     -  z  cos  0   ) 
3j   e     ^  ^  .  .  (5U) 

for  h  £  z  £   z      ' 

In  that  part  of  the  s       layer  betv/een  the   dipole  and  the  lov;er  "boundary 
surface  of  the  layer,    i.e.,    z,^£  z  <  h,    the  primary  radiation  is  absent  by 
hypothesis  and  we   are  left  v/ith   (b)  and  (c)  giving 

ik   (x  sin  9     cos  Cp    +  y  sin  0     sinCp    +  z  cos  9   ) 
g  =     Ai     e      ^  + 

ik  (x  sin  0     cos(ff  +  y  sin  0     sin  ^'    -  z   cos  9   ) 
Bi      e      ^  ^  .  (55) 

for   z     T    <   z   <  h. 
s-1   ~ 

Expressions    (5U)   and   (55)  give  us   the   form  of   our  plane  v/ave   function  g 
everywhere   in  the   s        layer. 

The  J        layer   (j   7^   s) 

In  any  layer  other   than  the   s      ,    g  evidently  consists   of  tv/o  plane  v/aves: 

(a)  A  wave   of  -Oiil-noun  amplitude  A^    traveling  with   the  primary   radiation. 

(b)  A  wave  of  Uiikno\-n  amplitude  BJ    traveling  against   the  primary  radiation. 
Thus  v;e  have 

.  .      ik.(x  sin  9"^    cos  (/"^   +  y   sin  6"^    sin  a*^   +   z  cos  0*^ ) 

g-^   =  J^   e     "^  '  + 


3,    e 


ik.(x  sin  0*^   cos  U."^   +  y   sin  9*^    sinj/;'^   -   z   cos   9*^) 
1  ^  (56) 

(j   =  0,1, 2,. ..s-1,    s+l..,.n+l) 


for  z .    ,    i  z  £  z .    . 
J-1  J 


37. 

(It  should  be  noted  that  A,  =  B,    =0  since  above  the  multilayer 

there  will  "be  no  plane  \/ave  moving  against  the  primary  radiation,  and  below  the 
multilayer,  i.e.,  inside  the  earth,  there  vrill  be  no  plane  uave  moving  \.'ith  the 
primary  radiation, )  .  n  - 

V/e  have  no\!   to  determine  the  unknovm  pnplitudes  A-.  and  Blr  through  the 

imposition  of  the  boundary  conditions  (3I)  and  (32)  of  Section  II.  These,  \;e  re- 
call are 

k^'  g^  =  i-.l'  g^*'  .  (57) 

^^=^Ldl    ,  (58) 

(j  =  0,1,2,  .  .  .  n) 

Since   conditions    (57)   ^^'^   (5^)   are   to  hold  at    the  surfaces   of   every 
layer,    v/e  now  allo\/  j   xo   range    from  0   to  n  including   the  values   s   -  1   and  s.      If  \fe 
apply   (57)   a^^d.   (5s)   to   (5^),    (55).    SJ^O-  (56)   and  observe   that,    as   in  Section  II, 
Snell's   lav;  is   a  consequence   of  the  boundary  conditions,   we   get  the   follov/ing  equa- 
tions for  A;      and  3^   : 

ik.z.   cos  9  .     -ik.z.    cos  0 

k^(b    .   +  A,^)   e      J   "^  +  k^  B^   e        ^   '^ 

J        so  1  J      1 

i+1  T+1 

„  .^,    ik...  z.    cos  0"  „  ._^,      -ik.^^z      cos  ©" 


ik.z.   cos  0  .      .     -ik.z.    cos  9 

k.    cos  ©*'(  8    .   +  A!')   e     "^   "^  -  k.   cos   9''   5^   e        "^   '^  - 

J  sj  1  J  1 

1+1  i+1 

■^i       -^T        xk..z.   cos  Q  .^,     .^-    -ik  ._^,  z  .cos  0*^ 

k.^,    cos  ©J"^  Af  ^      e     -^"^  J  -  k.   ,    cos   9"^"^"^        J^^  ^ 

J+1  1  J+1  1 

(60) 

(.j   =  0,1,2, ...n) 

0  ,  J  ^  s  0    n+1 

where  b  .  =  ,    •      and  where  A,  =  E,   =  0  . 
SJ    1  ,  J  =  s  11 


38. 

Equations  (59)  and  (60)  constitute  2n+2  equations  for  the  2n+2  uniinovm 
A^^  and  B^  ,  rememberinei'  that  all  the  Q*^s  can  be  expressed  in  terms  of  any  one  of 

them  (e.g., 9^)  throut^h  Snell's  lav/  (35)' 

'.vith  A''  and  B'^  deternined  by  (59)  ?r.d  (60)  v;e  may  regard  the  plane 
T,/ave  function  g  associated  "ith  the  first  subsidiary  problem  as  known,  and  \re   now 
solve  the  first  subsidiary  problem  by  the  use  of  the  superposition  method  of  Sectionll, 

The  solution  '/ill  be  given  by  (U5).   Since  o-x  dipole  is  nov/  at  the 

point  (0,0, h)  in  the  s   layer,  the  variables  of  integration  \rill   be  0  and  (p 

instead  of  Q  and  gj      (of  com'se,^  =  p   ).   lioreover,  the  function  A(0  )  given 

by  (1+9)  becomes 

ik    -ik  h  cos  G 
A(9^)  =  ^  e   ^  sin  0^  .  (6I) 

From  (U5),  (61),  and  the  expressions  (5^),(55)i  and  (56)  for  g,  the 
solution  of  the  first  subsidiary  problem  is  therefore  (AT  and  3^^  being  determined  by 
59  and  60)  "^ 


ik  -ik  h  COS0  ik  (x  sin  ©  cos  (p   +y  sin  0^sin  (/^^  +z  cos  0^) 

yWd^y-^   '      sin9'(l+A^)e  '  ^  d0^ 

0   0 
TT 

2n  2 

/  ik   -ik  h  cos  9  ik  (x  sin  6  cos(^^+  y  sin  0  sinc^  -  z  cos  Q^ ) 

d^V  -^  e   ^       sin  9^B^  e  ^  ^  ^  d0% 

0   0  (62) 

for  h  £  z  <  z 

s  • 

T  The  integration  contour  in  (62)  and  (63)  is  the  same  as  the  one  associated  \/ith  the 
integral  representation  of  the  primary  radiation.  This  follows  from  the  superposi- 
tion principle  as  expressed  by  (UU)  and  (U5)  of  Section  11.   Since,  in  the  first 
subsidiary  problem,  only  the  upper  half  of  the  dipole  radiates,  the  appropriate 
contour  is  that  of  (Uj).  Sae  Fig.  U. 


I .!      ''' 


-U 


39. 


ik.(x  sin9''cos(y5'+y  sin  0'^sin(/'+z  cos  G*^ ) 
L|;^=  /d^F/    -^e        "  sin  O''  A^le     J  ^  ^^  dO^ 

0       0 

2n       —  -1  00  .  .  .  . 

.,    ,  ^s  .    ik.(x  sin  Q'^cos&f^  +y  sin  0'^sin(>,'^   -z   cos  0"^) 

,     Ik     -ik  h  cos  9      .      qS_o       j  r       "  r  '   .s 

(63) 
0     u 

for  all   z   not    in   the  ran.re  h  1   z  -   z      « 

s 

Expressions    (62)   and  (63)   constitute   the   solution  of  our  first   sub- 
sidiary protlem,    remenibering   that  A,    =  3,        =  0, 


17.    Solution  of   the   Second  Subsidiary  Problem. 

Me   recall  that    in   the   second  subsidiary  problem   only   the  lo\;er  half  of  the 
dipole  radiates,    i.e.,   we   suppress   the  primary  radiation   in   the  region  h   <   z   <   z      , 
We  seek  the  ultimate   disturbance  U/     \inder   these   conditions. 

The  procedure   is    exactly  analogous   to   that  used  for  the   solution  of   the 
first    subsidiary  problem.     We   find  the   form  of  the  plane  \/ave  function  g   in  any   of 
the   layers   in  Fig.    5*      If  "e   still  associate  the   A-amplitudes  v/ith  plane  v;aves  mov- 
ing  in  the   same   sense  as   does   the  primary  radiation  and  the  B-amplitudes   with  plane 
waves  moving  against   the  primary  radiation,    we   see   tha.t   g  must  have  the   form 

ik  (x   sin  9^    cos  o/  +  y   sin  9^    sin  (/   +  z   cos  9^) 
g     =   (1   -t-  Ag)    e  +  .  ■ 

ik   (x  sin  9      cos  (fi     +  y  sin  9"    sin  (f    ~  z  cos  9    ) 

B2  ^  '  .        (6U) 

for  z  ,  <  z  <  h 
s-1  ~ 

.    .   ik.(x  sin  G"*  cos  (f'^+y   sin  9"^  sin  (fP  +  2.   cos  9"^) 

g-^   =  A^   e      J  ^  + 


ik.(x  sin  9*^   cos  ^   +  y  sin  9"^    sin  (/J^  -  z  cos  9*^) 

(j   =  0,1,2. ,n+l) 


^i^     ^  '  .  -  (65) 


for  all   z   not    in  the   range   z     -    *^   z   <  h 

s-1  ~ 


1+0. 

We  note  no\/  that    in   (65)  A^       =  B     =  0,    since  atove  the  multilayer 
there  will  be  no  plane  wave  moving  vath  the  primary  radiation,  \/hile  "belo\;  the 
multilayer,    i.e.,    inside  the   eartn,    there  will  "be  no  plane  \;ave  moving  against  the 
primary  radiation. 

Applying  the  boundary  conditions    (57)  and  (58)   to   (6U)   and  (65)  for  the 
determination  of  the  unknovm  amplitudes, we  evidently  require   (Snell'^s  law  is  again 
a  consequence  of   the  "boundary   conditions) 

.      ik  z .   cos  9"^  „      .     -ik.z.   cos  9"^ 

k^  A^  e     J   J  +  k^  B^  e     J   -^ 


i+1  i+1 

,   ,        ik.^,  z.cos   9^^  ^        .^^      -ik.^,  z.   cos  9*^ 


ik.z. cos  9*^  .    .     -ik.z.    cos  9'' 

k.    cos  9*^  A^  e     "^   "^  -  k.    cos  9-^3^  e       "^   "^ 


Vl°°=«'''<^3-l.j<'' 


.^,      ik._^,z.cos  9-^  .^T     .^,  ^k.^-z.    cos  0'' 

k       cos  ©Be  ^g^ 


(j   =  0,1, 2, ...n) 


°-^  ^    ^-^  ,n^l  0 


v/here    S,   1    ^  =  and  A^       =  3^     =  0  . 

^-^•J  l..i  =   s-1  2  2 


,j  =   s-1 

Once   again,   v;e  have   2n+2  equations    (66)   and  (67)   for   the  2n-*^  unknowns  Ap  and  T^^   , 

1  s 

remembering  that   all   the  9  s  can  be  expressed  in  terms  of  9     through  Snell's  law  (35). 

V/ith  the  amplitudes   thus  determined,   \/e  may  write  dovm  the   solution  of 
the  second  subsidiary  problem.      V/e  again  use   (U5)  as   our  basic  expression.      Since 
vie  are  now  considering  that   only  the  lov/er  half  of  the  dipole  radiates,    however,    the 
integration  contour  appropriate  to  the  solution  v/ill  be  that   of  expression  (U8)  in 
Section  II.     This    contour  is  also  shovm  in  Fig.    U. 

The  function  A(e    )   being  given  by   (6I),   the  function  g  by  (6U)  and  (65), 
and  the  amplitudes  A^  and  Bp  being  determined  by  (66)  and   (67),  we  have   as   the 
solution  of  the  second  subsidiary  problem- 


.-s,-...  -'Mi:  .'  ■■  !' 


1 1«       •■   '  -• 


•Tii''  "I 


Ul. 


2n      TT     . 

fik     -ik  h  cos   G^  ik   (x   sin  0^cosaf+  y   sin  O^sinr^-    +  z  cos  0") 


sin  e   (1+A^)e 


60" 


2n    TT 


ik     -ik  h  cos  0 


+  ioo 


ik  (x  sin  0  cos^+  y  sin  9  sin/p  -  z   cos  0    ) 


dO^, 


(68) 


for  z     -    <   z   <  h 
s-1 


ik.(x  sin  0'^cosa-'^+  y  sin  e'^sin65'+  z  cos  Q"^) 


d©' 


2n    V 


ik       -ik  h  cos   0 
^/d<^/  "2^  e        ^  sin  e^B^  e     -> 


0       ^TT   ^    . 

p    +    ICD 


ik.Cx  sin  ©■^cosa>^+  y  sin  9'^sinrt?''-  z   cos  0"^)        ^ 


(j  =  0,1,2, ...n+l) 


(69) 


for  al  z  not  in  the  range  z  ,  £  z  <  h  . 

Sxpressions  (68)  and  (69)  constitute  the  solution  of  the  second  subsidiary 
problem,  remembering  that  Ap  =  Bp  =  0. 

18.  Conclusion  of  Section  III. 

The  solution  of  the  propagation  problem  of  Fig.  5,  i.e.,  the  problan  of 
the  dipole  within  a  multilayer  region,  is  given  by  ^J^   +  \^ ^,   \heTe  \\f-^   and  \^^   are 

the  solutions  of  the  first  aiid  second  subsidiary  problems,  respectively.   If  we 
carry  out  the  addition  of  K\)-.    and  \0      it  is  easy  to  see  that  the  final  solution^  of 
the  TDroblem  of  Fig.  5  "i^^  ^^ 


U2. 


^^s^      /       /ik       -ik  h  cos  0^  .    ik.(x  sin  0'^cos6&^+  y  sin  Q'^sinaP+z  cos  &) 

/      p.  I       o.  <^  .5.11  T         "  T 


V=fes,^R-V**'M'      =  slne^Aje    J  ^  ^  «= 


/     ik       -ik  h  cos  0^  ,      ik.(x  sin  9"^cos(V^+  y  sin  Q'^sinc^-z  cos  8^) 

d^y     -2^  e        ^  sin  Q%{   e     ^  ^  ^  d©^ 

(70) 

ik       -ik  h  cos  9^  ,    ik.Cx  sin  Q'^cos^+  y  sin  9'^sinQp+2  cos  9"^) 

sm  9  Ape     *^  d9 


ik       "ik  h  cos   ©^  .    ik.(x   sin  9"^cosar+y   sin  9'^sin^?'-   z   cos  ©"^ ) 

+  /d^/_ie       ^  sin  9^B^e     J  ^  '  d©^ 

(J  =  0,1,2, ...n+1)   ■ 


0.  j  f 

1,  J  =  s 


\?here  ^  •  =  .A?  and  B''  satisfy  (59)  and  (6o), 


A^  and  3^  satisfy  (66)  and  (67),  aJ  =  ^^"^     =  k^^=   3°  =  0,  </>.  =  (^^  for  all  j,  and 

and  R  is  the  distance  from  the  observation  -point  to  the  point  (0,0, h)  \/here  the 

•t  _  _     _ 

dipole  is  located,  xhe  solution  vectors  S  and  H  maj-'  now  he  derived  from  Uj  through 

relations  (5)  and  (6)  or  (S)  and  (9)  of  Section  I. 

It  should  he  noted  that  the  usefulness  of  the  sclation  (7O)  depends 

entirely  on  ho\/  successfully  one  evr.luates  the  complex  integrals. 


'pAccording  to  the  rule  explained  in  Section  II,  in  the  first  tuo  integrals  of  ex- 
pression (70)  0  <  Re  e*^  <  ^  while  in  the  last  tv/o  integrals,  H  <:  Re  Q J  <  n  . 


^3. 

Section  IV 
A  Special  Case  -  n  =  1. 

19 .  Statement  of  Proplem. 

In  this  section  v/e  shall  obtain  the  solution  of  a  special  case  of  the 
propagation  protlem  of  Fig.  5  i>^-  Section  III.   The  special  case  is  the  one  for  which 
n  equals  1,  i.e.,  referring  to  Fig.  6,  our  problem  is  as  follovis: 


/^  5c€,  CO  -i-lciUJ 


G 


R 


t 


Fig.  6 

The  atmosphere  is  divided  into  tv;o  layers,  the  first  of  height  z,  with 
constants  £  and  C"  ,  and  the  second  of  infinite  height  with  constants  £   and  <X>. 

A  vertical  electric  dipole  A  of  ar^ilar  frequency  ^  is  situated  v/ithin  the  first 
layer  at  a  height  h  above  the  earth.  We  seek  the  solution  field  ^^ anyi/here  on  or 
above  the  earth. 

To  solve  this  ;jro":;lem  v.e   make  direct  substitutions  into  the  equations 
and  formulas  of  Section  III.  Specifically,  we  solve  (59)  and  (60)  for  the  amplitudes 
associated  v/ith  the  first  subsidiary  problem, (66)  and  (67)  for  the  amplitudes  associ- 
ated v:ith  the  second  subsidiary  problem,  and  substitute  these  results  inT;o  (70)  to 
obtain  the  final  solution. 

We  now  carry  out  these  operations. 


20.  Solution  of  the  First  Subsidiary  Problem. 

Our  equations  are  (59)  a^^cL  (60)  specialized  for  n  =  1.  Remembering 

0         2 

that  A-    =  B,    =0,    these  become 

11  2 

ik,  z   cos  9  „  ,      -ik  z  cos  G  ^^o^     cos  © 

kj(l  +  Aj)e     ^   ^  +  k^J  e        ^  ^'  =  k|  Aj  e     2   1  ^  (yi) 

1  1  2 

ik  z  cos  9  ,  -ik^z   cos  9  ik  z   cos  9 

k^   cos  e   (1  +  Aj^)e  -  k^cos  9  B^   e  =  k  cos  9'^A^  e     ^ 

(72) 

2     0  2     1  2     1 

kjB^     =k^  AJ:  +  k^B^      ,  (73) 

-k^  cos  9°  B°  =  kj^  cos  9''"  A^  -  k^  cos  9''"  B?"  (yU) 

The  exponentials   are  absent   from  equations    (73)   and  (7U),   because   in  equations    (59) 
and  (60),    Zq  =  0, 

Equations   (7I)   through  (7U)  constitute  four  equations  for   the  four  un- 
knovms  A^    ,  A^    ,   B^^  and  B,    .   To   solve  them,   ve  proceed  as    follov/s:     we  multiply 

equation  (7I)  by  cos  9     and  equation   (72)  by  k     and  subtract  equation   (72)   from 

equation  (7I),    resulting  in 

12  2  1     ik  z  cos  9  -ik  z  cos  9-'' 

Aj^(k^  cos  9     -  k^k^  cos  9   )e         "•  +  B^(k^  cos  9Sk^k     cos  9^)   e 


,         P  P     ik^  z^  cos  9 

(k^k^  cos  9  -  k^cos  ©    )e 


(75) 
We  no\/  multiply  equation   (73)  by  cos  9     and  equation   (7U)  by  k^,   add  and  obtain 

A^(k^  cos  ©     +  k^k^cos   9-'-)  +  B^{k^  cos  9°  -  k^k^  cos  9-'-)     =     0.  (76) 

Cancellation  of  a  k,  from  both  sides  of  equations  (75)  and  (76)  results  in 

1,       ?       1  ik.z^cos  9  -ik^z  cos  ©■'• 

A^(k^cos  ©  -k^cos  9  )e    "•      +  Bj^(k^cos  9'^+  k  cos  9  )e 


ik,  z,  cos  9 


1 


1  2     ^h  . 

(k^  cos  9  -  k^cos  9'^)e  (77) 

A^(kj_  cos  9*^  +  k^  cos  ©■'■)  +  BJ^(k^  cos  9°  -  k^  cos  9^)  =  0  .  (78) 


I/.   ;. 0 


Equations  (77)  ^^^   (78)  can  now  "be  solved  for  A,  and  B,  .   Solution  of  equation  (78) 

for  A^  gives 

k  cos  9  -  k^  cos  9 

h   "   k^cos  91  +  k^  cos  9U   ^1  =  ^0^1   •  ^79) 


uhere  f^  = 


1   ,      ^0 

0 


k^  cos  9  -  k^  cos  9 


0    ,      .1  .  ,      .0 
0 


k^  cos  0  +  k.  cos  9 


Substituting  from  equation  (79)  into  equfition  (77)i 

,        P  .   ^^l^^l  cos  9  -   -ik  z  cos  9 

Tq  B^(k^  cos  9  -  k^  cos  9  )e  +  B^(k^  cos  9  +  k  cos  9  )e 

ik, z, cos  9 


1  p   -'^i ''-.^^^  ^ 

(k  cos  9  -  k^  cos  9  )  e  ,   (80) 


iifhence  v/e  obtain 


,  P   ik  z  cos  9 

,  (k  cos  9  -  k  cos  9  )  e 

B^=  . ^ i ^ 

2  1   i^z,cos  9  -ik  z-cos  9 

f  (k^  cos  9   -  k  cos  9  )  e  +  (k.^cos  9  +  k  cos  9  )e 

^         ^  (81) 


1  2 

k  cos  9  -  k.  cos  9 

V/riting  f^  =  y —   ,  v;e  finally  get 

k  cos  9  "*■  k^  cos  9 


r.        2ikT  Zt  cos  9 

^1      ^1  ^   ^  ^ ,   , 

1  ~  2ik  z  cos  e^     •  ^^^•' 

1-f,  f.e 

From  equations  (79)  and  (82)  \;e  obtain 

2ik^  z  cos  9 

a1=   ^Q  ^1  ' ,  (83) 

1  2itu  z,  cos  e       * 

^  -  ^0^1  ^ 


U6„ 


and  from  equations   (71),    (82)   and  (S3)   there  results 


2  ^1   2 

^1  =    ^4^ 


1  + 


Vi^ 


2ilc^  z      cos  Q 


1   -  f^f^   e 


2ik  z  cos~Qi 


i(k  z,    cos  Q 
.       1  1 


kpZ,    cos  0    ) 


k-i  f,    e  -i(kTZ^cos   Q  +k^z-COsGr) 

*  ^iT^  — ^^ r  ' 

2  2ik^z^cose  (g^j 

^0   1   ® 

Eouations  (82), (83),  and  (8U)  constitute  the  solution  of  the  first  sub- 
sidiary problem.  We  shall  not  bother  to  solve  for  B^ ,  since  we  are  not  especially 
interested  in  the   field  uithin  the  earth, 

21.    Solution  of  the  Second  Subsidiary  problem. 

In  this  case  our   equations   are   (6b)  and  (67)   specialized  for   n  =  1, 

2         0 
Remembering  that  Ap  =  Bp  =  0,    these  become 


2     1      i^^z^cos  9  "^^1^1   ^°^  ®  ?     2     ~'^^?^1  ^°^   ®' 


k^  A^  e 


+  k^  B^  e 


k^B^e 


(S5) 


-    ik  z  cos  0       •        -ik  z  cos  © 
k^  cos  9^k^  e         ^  -  k^  cos  6  3^  e   ^  ^ 


k2  A°  =  k^  (1  ^  A^)  *  k^  B^  , 


P  p  -ik  z  cos  © 
■k  cos  Q'^'Z   e   ^  -^ 
^      ^  (86) 


(87) 


k^  cos  0  A°  =  k^  cos  e"""  (1  +  A^)  -  k^  cos  ©'^  Bp 


0 


izs) 


0  .1 


Equations  (85)  through  (88)  constitute  four  equations  for  the  foxir  unknovms  A  ,  A  , 


Bp  and  3p  , 


2'   2' 


H7. 

2 
I'ialtipl^ine   equation   (85)  "by  cos  6     and  equation   (S6)   by  k     and  adding, 

we  obtain 

1       ilS  V°^  ®       1     2  2  1     -i^2  ros  6 

A^CkJ  cos  Q'^  +  k^k^  cos  9   )   e  +  B^lk^  cos  6''  -  k^k^cos  0   )e  =0. 

(89) 

Multiplication  of  equation  (87)  by  cos  0  and  of  equation  (88)  by  k  and  subtraction 
of  equation  (88)  from  equation  (87)  results  in 

A^(k^  cos  9°-k  k  cos  ©■"")  +  ^\^^'l   cos  9°  +  k^k^  cos  9''")  =  k^k^  cos  9^  -  k^  cos  9^.    ; 

(90) 

Cancellation  of  k,  from  both  sides  of  equations  (89)  and  (90)  gives 

,        P  ,    ik^z^cos  9    ^        P  -ik.z-cos  9 

Agdc^  cos  9  +  k^  cos  9  )  e  +  ^2^  '^°^         "  ^2   ^°^  ^^®  ""  °* 

(91) 


^(k^  cos  9     -  k^cos  9^)  +  B^ik^  cos  9°  +  k^  cos  9"'')  =  k^  cos  ©"^  -  Ic.    cos  9°       (92) 


A2 

Equations  (31)  and  (92)  can  nov;  be  solved  for  Ap  and  3  , 
Solution  of  equation  (91)  for  Ap  gives 

12  1  1 

^  k     cos  9  -  k     cos  9         -2ik^z,cos  9     ,  -2ik^  z     cos  9 

4  = r-^ r  ^  ^2     =h^  ^2'       ^53) 

k     cos  9     +  k^cos  9^^  ^-  x  d 

Substitution  of  expression  (93)   into  equation  (92)  gives 

1  n  1        -2ik^z^cos  ®  1  r,  T 

B^  f^    (k^  cos  9     -  k^cos  9   )   e  +  B2(k^cos  9  H-k^cos  9   )  =  k  cos  e^k^cos  9° 

(9U) 
whence  we  obtain 

,                        k^  cos  9.  -  k,  cos  9 
bJ  =  '^ 1 1 .     (95) 


-2ik^z,  cos  9 


f^(k^  cos  9°  -  k^  cos  e""-)  e      ^       +  k^  cos  9*^+.  k  cos  9'' 


1(S. 


B^  =■ ^2 _  (56, 

-2ik^  z  cos  9 
1  -  f^f,  e 


From  equations  (93)  ^'^^   (96)  there  results 


-2ik^  z^  cos  9 
,1       ^0^1  " 


^2  =  r; TT-  (97) 


-2ik.  z  cos  © 
1  -  f^f^^  e    ' 


and  finally  from  equations  (85),  (93)  and  (96)  we   have 


2        "^^^'-1^1  "^^  ^  2  1 

K.  i„f,  e  i(k_z,  cos  0  +  k,  z.  cos  ©  ) 

B'^  =   (— )  ^L± e   ^  •'•  -'■•'■ 

?      k  1 

2  ~2ik  z   cos  0 

1  -  f,f^  e 


2  2         1 

k,  f_  i(k„ZT  cos  9  -k,  z^cos  ©  ) 

+  (J.)     . — 2 e   ^  ■'•        ■'•■'■ 

k  1 

2  -2ikT  z-  cos  ©  ,^_. 

1  -  f  f  e    ^  ^  (98) 

0  1  ^ 

Equations  (96),  (97>  and  (98)  constitute  the  solution  of  the  second  subsidiary 
problem. 

With  the  first  and  second  subsidiary  problems  solved  for  the  special 
case  of  the  propagation  problem  shown  in  Fig.  6,  v;e  may  v/rite  do\/n  the  solution 
of  this  problem  in  accordance  uith  expression  (JO)  of  Section  III.  This  solution  is 
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